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The ^Sq, '^Si, and ^Di nucleon-nucleon scattering phase shifts are cal- 
culated at next-to-next-to-leading order (NNLO) in an effective field theory. 
O I Predictions for the ^Pi, ^^0,1,2, ^D2, and ^-D2,3 phase shifts at this order are 

also compared with data. The calculations treat pions perturbatively and in- 
clude the NNLO contributions from order and Qf radiation pion graphs. 
In the ^/Si, '^-Di, and ^-Po,2 channels we find large disagreement with the Ni- 
jmegen partial wave analysis at NNLO. These spin triplet channels have large 
I corrections from graphs with two potential pion exchange which do not vanish 

in the chiral limit. We compare our results to calculations within the Wein- 
berg approach, and find that in some spin triplet channels the summation 
of potential pion diagrams seems to be necessary to reproduce the observed 
phase shifts. In the spin singlet channels the nonperturbative treatment of po- 
tential pions does not afford a significant improvement over the perturbative 
approach. 
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I. INTRODUCTION 



Understanding how nuclear forces emerge from the fundamental theory of Quantum 
Chromodynamics (QCD) remains an outstanding problem in theoretical physics. To study 
the physics of hadrons at scales where QCD is strongly coupled, it is useful to employ the 
technique of effective field theories. Model independent predictions for low energy nuclear 
phenomena can be made by using an effective Lagrangian which includes nucleons and 
pions as explicit degrees of freedom and all possible interactions that are consistent with 
the symmetries of the underlying QCD theory. This method, known as chiral perturbation 
theory, has been successfully applied to processes involving and 1 nucleons (see e.g. 

Weinberg |Q originally proposed using effective field theory for few body problems in 
nuclear physics. Weinberg's procedure applies ordinary chiral perturbation theory power 
counting to the nucleon-nucleon potential and then solves the Schrodinger equation using this 
potential. Phenomenological studies of NN scattering phase shifts and deuteron properties 
which use this technique can be found in Refs. 

Application of effective field theory to two nucleon systems is complicated by the exis- 
tence of a shallow bound state in the spin triplet channel and the large scattering length 
in the spin singlet channel. In Refs. Kaplan, Savage, and Wise (KSW) proposed a 

new power counting which accounts for these effects. This approach is more like ordinary 
chiral perturbation theory in that power counting is applied to the amplitude rather than 
the potential. All observables are expanded in powers of Q/A, where Q is either or 
p (the nucleon momentum), and A is the range of the effective field theory. Because the 
S-wave scattering lengths (denoted by a) are large, powers of pa must be summed to all 



orders |jTO[|. This requires a nonperturbative treatment of the leading 4 nucleon operators 
with no derivatives. Higher derivative operators and pions are treated perturbatively. The 
perturbative treatment of pions makes it possible to obtain analytic expressions for ampli- 
tudes. One theoretically appealing aspect of the KSW power counting is that all ultraviolet 
divergences appearing in loop graphs are cancelled by contact operators appearing at either 
the same or a lower order in the expansion. This is in contrast with Weinberg's approach 
in which unsubtracted divergences introduce cutoff dependence which is cancelled at higher 
orders in the expansion. The residual dependence on the cutoff gives an estimate of the size 
of higher order corrections. 

It is clear from naive dimensional analysis that the KSW expansion will converge slowly. 



2 



To see this, compare the contribution to the amphtude from single pion exchange and the 
pion box diagram: 



J— 2/^ \m, 



where A, B are dimensionless functions. The factor of M, the nucleon mass, comes from 
performing the energy integral by contour integration and taking a pole from one of the 
nucleon propagators. The factor of l/(47r) is an estimate of the size of the loop correction. 
(If a pion pole is taken the contribution is smaller by m^/M From Eq. (P one expects 

an expansion parameter of order {g\m.„M) / (Sir p) = ni-^/ANN — 0.5 0. This suggests that 
perturbative pions will converge, albeit slowly. 

Many processes involving two nucleons have been computed to next-to-leading order 



(NLO) in the KSW expansion |12,13|. The results of some of these calculations are reviewed 
in [|1^]. Typically, one finds 30%-40% errors at leading order (LO) and 10% errors at 
NLO. These results suggest an expansion parameter Q/A ~ 1/3 or A ^ 400 Me V. This 
is consistent with the estimate of the expansion parameter given above. Obviously, it is 
important to extend existing calculations to higher orders to see if the convergence of the 
expansion persists. 

At the present time, few NNLO calculations^ are available in the theory with pions. The 



deuteron quadrupole moment is calculated to NNLO in Ref. [T^. The result of the NNLO 
calculation of the ^5*0 phase shift has been presented in Ref. [0,|l3l and independently in Ref. 
18| . However, these ^Sq calculations are incomplete because the full order Q contributions 



from radiation pion graphs were not included. The ^Si — ^Di mixing parameter is calculated 
to NNLO in Ref. [|19|, where it is demonstrated that the expansion is converging for p < 
140 MeV. For these momenta, the error is comparable to that of calculations within the 
Weinberg approach 0. 

In this paper, we present NNLO calculations of the ^5*0, ^5*1, and ^Di phase shifts in 
nucleon-nucleon scattering, including contributions from radiation pions. At this order we 
find that the radiation pion diagrams have trivial momentum dependence and their effect 
cannot be distinguished from the contributions of a local operator. In the ^5*0 channel, 
the NNLO fit agrees with data to < 1% accuracy for p ~ m.„. In this channel, the KSW 



^For the NN scattering amplitude LO is Q ^, NLO is Q'^, and NNLO is order Q. In this 
paper this terminology will be used even for cases where the LO contribution vanishes. 
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expansion works as expected. However, in the spin triplet channel we find that the expansion 
breaks down at NNLO. For the ^Si and ^Di phase shifts, the NNLO calculation actually 
does worse at fitting the data than the NLO prediction. In the ^Si channel, the NNLO 
corrections are as large as the NLO corrections for p = itlt^. In the ^Di channel there is no 
sign of convergence for any value of p. We find that the failure of the EFT expansion in these 
two triplet channels is due to large non-analytic corrections that grow with p coming from 
graphs with two potential pions. These terms do not appear in the spin singlet channel. The 
reason for the difference in the quality of the perturbative expansions in the two channels 
is that the potential between nucleons arising from pion exchange is much more singular in 
the spin triplet channel than in the singlet channel. We elaborate on this point in section |^ 
of the paper. 

Next, we examine the NNLO predictions for the P and D wave phase shifts drawing on 



results from Ref. [|rT|. At LO these phase shifts vanish. In these channels the only contri- 
butions at NLO and NNLO come from potential pion exchange. Contact interaction and 
radiative pion contributions do not enter until higher order. Thus predictions for these phase 
shifts contain no free parameters, and it is possible to unambiguously test the perturbative 
treatment of pions. In the spin singlet channels (^-Pi, ^-02) corrections from two potential 
pion exchange are small, and the errors at p = are (13%, 33%). At p = m,r; the NNLO 
predictions for the ^Pi, channels have errors of the expected size (15%, 8%). In the ^Pq, 
^P2 channels errors are bigger than expected (170%, 52%). Like the ^5*1 and ^Di channels, 
these spin triplet channels have non-analytic contributions that grow with p. 

Our final section includes a comparison of our calculations with those of Refs. ||^,|^ which 
use the Weinberg approach. In spin singlet channels the corrections obtained by summing 
perturbative potential pion exchange to all orders are negligible. In particular, in the ^Pi and 
channels, single pion exchange gives the same answer as the LO Weinberg calculation 
which treats potential pions nonperturbatively. Here corrections from soft and radiation 
pion graphs as well as contact interactions appear to be much more important. In the KSW 
expansion, these effects appear at one higher order than the results presented in this paper. 

In some spin triplet channels (^S*!, ^Pi, ^-Di) the summation of potential pions gives 
significant improvement relative to the calculation which treats the pion perturbatively. 
There are also spin triplet channels where nonperturbative potential pions seem to be less 
important than soft pion graphs and four nucleon operators. This is true in the ^Pq and 
^P2 channels, where the LO calculation in the Weinberg scheme does no better than the 
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LO term in the KSW expansion. Finally, in the ^-D2,3 channels, the KSW expansion at 
NNLO gives predictions that are as accurate as the NNLO Weinberg calculations and so a 
nonperturbative treatment of pions does not seem to be necessary in these channels. 

The rest of the paper will be organized as follows. In section II, the formalism relevant 
for our calculation is introduced. We define all operators appearing in the Lagrangian to 
the order we are working and discuss the solution to their renormalization group equations 
(RGE). Solving the RGE perturbatively ensures that observables are renormalization scale 
independent, as in pion chiral perturbation theory. We also discuss our method for fitting 
the constants at each order in the expansion. In section III, expressions are presented for 
the ^5*0, '^5'i and ^Di amplitudes up to NNLO. Detailed comparison of the theoretical phase 
shifts with the Nijmegen phase shift analysis pO|] appears in this section. In Section IV, 
we look at NLO and NNLO contributions to nucleon-nucleon scattering in the ^Pi, ^-Po,i,2 5 
^D2, and ^-D2,3 waves. In the final section, we discuss our results and their implications 
for the perturbative treatment of pions. Details of the calculations are contained in the 
Appendices. In Appendix ^ we describe a trace formalism for projecting partial wave 
amplitudes from Feynman diagrams. In Appendix ^ we give explicit expressions for all 
individual graphs at NNLO, except for graphs involving radiation pions. We also describe 
a general strategy for analytically evaluating massive non-relativistic multi-loop Feynman 
diagrams. In Appendix |C|, the S-wave radiation pion contribution is discussed in detail. 
The power counting for radiation and soft pions is reviewed and the complete order Q 
contribution is evaluated. 



II. FORMALISM 



In this paper, we will follow the notation in Refs. [p|, p!9| , pT| . The relevant Lagrangian for 
NN scattering at NNLO is 



C = ^Tr + ^ Tr(m,S + m,St) + A^t (^tDo + ^ ) ^ 



^9A 
2 

64 



+ f^^Tr(m«)0(^) - 



D 



(^) 

4 2 
UJ 



{Tr2(m«) + 2Tr[(m«)2]}oJ 



(2) 



5 



Here qa = 1-25 is the nucleon axial-vector coupling, S = = exp(2in//) where 



n 



TT 



TT 



v 



TT 



-TT 



(3) 



/ = 131 MeV is the pion decay constant, the chiral covariant derivative is = 5^ + 
+ i'^d^^), and = \{^mq^ + C.^rriq^'^), where rriq = diag(mu, m^) is the quark mass 
matrix. At the order we are working co'Tr(m^) = uj^rriu + rrid) = m"^ = (137 MeV)^. In 
Eq. (0), s = ^Sq or ^^i. Below this superscript will be dropped when it is clear from the 
context which channel is being referred to or when the reference is to both channels. The 
two-body nucleon operators are: 



o. 



(^) 



tT p(s) 



tT p(s) 



of ^) = (iV^p/'^^)iV)t(iV^pf ^^)iV) + h.c. 



(4) 



where the projection matrices are 



P> 



iy/n - 1 ^ n ' ■' 



(5) 



and V = The derivatives in Eqs. (1) and (|^) should really be chirally covariant, 

however, only the ordinary derivatives are needed for the calculations in this paper. 

Ultraviolet divergences are regulated using dimensional regularization. All spin and 
isospin traces are done in n dimensions, where c? = n -|- 1 is the space-time dimension. 
Regulating the theory in this way preserves the chiral and rotational symmetry of the theory 
as well as the Wigner symmetry [0,^ of the leading order Lagrangian, as discussed in Ref. 

Hi- 

The KSW power counting is manifest in renormalization schemes such as power diver- 
gence subtraction (PDS) | 



I or off-shell momentum subtraction (OS) p4| , p5| , pT| . (In this 
paper the PDS scheme will be used.) In these schemes the coefficients of the S-wave opera- 
tors in Eq. (|) scale as ~ l/(MA"/i"+i), where \i is the renormalization scale, and A is 
the range of the effective field theory. The renormalization scale is chosen to be on the order 
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of the nucleoli momentum p which is of order m,r- Letting jj, ^ p ^ m^r ~ Q the scahng of 
the coefficients in Eq. is: 

LO: Clf\f,)^l/Q (6) 
NLO: p'ci'\f,)^Q\ mlD^\^,)^Q' 
NNLO : / C^^ (/i) ~ g, mlp" (/x) ~ Q, ml D^^ (/i) ~ Q , p'cf (/x) ~ Q • 

Note that from simple dimensional analysis one would expect these coefficients to scale as 
C2n ~ l/(MA2"+i). However, these coefficients are larger than naive dimensional analysis 
predicts because the theory flows to a non-trivial flxed point for a —>■ ±00. (See Refs. P, |26|j27[| 
for a more detailed explanation.) Since Co'^^(/i) ~ l/Q, and each nucleon loop gives a factor 
of Q, power counting demands that graphs with Cq's be summed to all orders. This sums 
all powers of ap. Operators with derivatives or insertions of the quark mass matrix scale as 
Q^, n >0, and are treated perturbatively. 

In Eq. (0) we have not included four nucleon operators for partial waves with L > 1 
because these operators enter at order or higher. For example, the coefficients of the 
four P-wave operators with two derivatives are not enhanced by the renormalization group 
flow near the flxed point and therefore scale as 1/{MA^). Thus, these P-wave terms in the 
Lagrangian are order Q'^/{MA^). As a result the order Q predictions for partial waves with 
L > 1 come completely from pion exchange and have no free parameters. 

There is another term in the Lagrangian in Eq. with an S-wave four-derivative oper- 
ator distinct from 0'^\ C = C^^ O^^ where 



is) 



For the process N{pi)N{p2) —>■ N{p^)N{p4) this operator vanishes on-shell since energy- 
momentum conservation gives {pi — ^2)^ = {ps — PaY- In deriving the RGE's only on- 
shell amplitudes are relevant. In fact the off-shell Green's functions do not have to be ix 
independent, as illustrated by the off-shell C2 amplitude given in Eq. ( P28|) of Appendix ^ 



Thus, to derive an RGE for 6*4^ ■* (/i) it is necessary to consider an on-shell process in which 
this coefficient gives a non-zero contribution. Although C^^\^) does not contribute to NN 
scattering, it may contribute to interactions with photons when the operator in Eq. (0) is 
gauged. Diagrams with two 0^2'' operators renormalize O]^^ making C^^\^) ~ The 

is) (s) ~ is) 

fact that Oi rather than 0\ —0\ has an enhanced coefficient differs from the conclusion 



in Ref. 28 . 
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Relativistic corrections contribute at order Q to the S-wave amplitudes. They are sup- 
pressed relative to the leading order amplitude by {Q/MY rather than (Q/A)^. In Ref. 
p8| these corrections are computed and found to be neghgible relative to other order Q 
contributions. Therefore they are left out of our analysis. 

For momenta p ^ m.^ pions should be included in the theory. There are three types of 
contributions from pions: radiation, potential, and soft. In evaluating non-relativistic loop 
diagrams the energy integrals are performed using contour integration. When the residue 
of a nucleon pole is taken the pion propagators in the loop are potential pions. When the 
residue of a pion pole is taken the pion will be either radiation or soft. Potential pion 
exchange scales as Q", and is therefore perturbative. Radiation and soft pions begin to 
contribute at order Q and respectively. The power counting for pions is discussed in 
detail in Appendix 0. Because pion exchange is treated perturbatively the dominant scaling 
of the C2n(/i) coefficients is the same as in the theory without pions. 

The Q scaling in Eq. can be determined by computing the beta functions for the four 
nucleon couplings appearing in Eq. to the order we are working. The procedure used 
for computing beta functions in the PDS scheme is described briefly in Appendix B and 
in detail in Ref. Our results are slightly different than Ref. |^ because all spin and 
isospin traces are performed in n dimensions rather than 3 dimensions. For the ^5*0 channel, 
the beta functions to NNLO are: 



oil 
oil 



Oil 

oil 



Mil 
Arc 

-( 

Mil 



C0C2 

C0D2 1 1 



24^ + 3 



1 + 2 



2/2 An 
9a Mh' 



' 9a Mil' 
2/2 An ^ 



2/2 An J 
,9^Mii' 
'2/2 An I ' 2/ 



+ 



9\ 



My 
A^J 



(Co) 



An 
Mil 



877/2 

2C^C, + {C2f 



An 



An 
Mil 



2 Co-E'4 + 2D2 C2 + 2C2 Cq 



D2 ' 



9l 



An 



My 



2/2 \AnJ 



9l 



2C0D4 + W +2^2C7o|^(^^)" 



(8) 



where the contribution (3^^ from radiation pions is given in Eq. ( |(J38| ) . All coupling constants 
are functions of /i. In the ^5*1 channel, the beta functions for Co, C2, and D2 are: 
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and the beta functions for C4, E4 and are identical to those in the ^Sq channel. (We 
have corrected a sign error in the C2 ^^"^ beta function computed in Ref. [^].) The running 



of Cf^(/x) is discussed in Ref. [|T^]. Terms in the beta functions that vanish as /i are 
from linear power divergences and are renormalization scheme dependent. These terms are 
necessary for a consistent power counting near the a —* ±00 fixed points. Taking Qa ^ ^ 
gives the dominant power contributions, and these terms are the same in renormalization 
schemes with a manifest power counting like PDS or OS. Finally, terms that do not vanish 
as /i correspond to logarithmic divergences and are scheme independent. 

It is desirable that the amphtude, and hence all physical quantities, like the scattering 
length, be /i independent at each order in the expansion. This can be accomplished by 
expanding the coupling constants in Q |p5| : 



^0 ^ ^0 + ^0 -r ^0 

D2^D2 + D^i^^ . (10) 

The first piece of Cq is treated nonperturbatively (i.e. Cq ~ Q~^), while Cq^"* ~ Cq'^ ~ Q. 
Because of the perturbative expansion of the couplings in Eq. (|T0|) there are ten constants of 
integration that appear in the calculation of the NNLO S-wave phase shifts. However, the 
NNLO amplitude depends only on six independent linear combinations of these constants. 
The coupling constants are also subject to two further constraints: 

1. At this order, C4, i?4 and D4 are determined entirely in terms of lower order couplings 
as a consequence of solving the RGE's and applying the KSW power counting. 

2. Spurious double and triple poles in the NLO and NNLO amplitudes must be cancelled 
in order to obtain a good fit at low momentum. 
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X 



FIG. 1. Order 1/Q (first row) and (second row) diagrams for NN scattering. The solid lines 
are nucleons and the dashed line is a potential pion. 

An example of constraint 1 is provided by the solution of the RGE for C4 given in 
Eq. (D i: 



C4 



M 



(11) 



Co ' '"47r 

where p is a constant of integration. In the theory without pions, p is proportional to the 
shape parameter, which is ~ in the KSW power counting. In the theory with pions 
p ^ too, since its size is determined by the scale A. Therefore, {C2Y /Cq ~ Q~^, while 
p(Co)^M/(47r) ~ . The second term is subleading in the Q expansion, and should be 
omitted at NNLO, so C4 = {C2Y /Cq. Solving Eq. (^ gives similar relations for £'4, and D^. 



2C2D2 



+ 0{Q- 



{D2 



0{Q 



-2\ 



(12) 



assuming that the constants of integration are order The beta functions for and D4 
have contributions from chiral logarithms, which are determined by the ln(/i) in D2- 

Constraint 2 is due to the nonperturbative treatment of Co, which gives rise to spurious 
poles at higher orders in the expansion. The leading order amplitude A-i has a simple pole 
at p = ^7. The NLO amplitude is proportional to J^ii, and therefore has a double pole, 
while the NNLO amplitude has terms proportional to and J^_i. To obtain a good fit at 
low momentum, parameters need to be fixed so that the amplitude has only a simple pole 
at each order in the expansion. This requires that A-i have its pole in the correct location 
and that the residues of the spurious double and triple poles vanish. This requirement leads 
to the following good fit conditions [p5|]: 



1 



^-1 







^0 



0, 



^1 



p=p'' 



0, 



(13) 



where p* is the location of the pole. The second condition first appears at NLO, the third at 
NNLO. The residue of the triple pole in Ai vanishes by the second equation in Eq. (|l^). The 
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FIG. 2. Order Q contact interaction and potential pion graphs for the and ^Si channels. 
At this order the first three graphs do not introduce new parameters as explained in the text. 
Radiation pion diagrams with order Q contributions are shown in Appendix 

first equation results in 7 = —ip*, while the other equations give constraints which eliminate 
two of the remaining parameters. In order to solve the constraints in Eq. ( |TB| ) we must allow 
the coupling constants Cq*^ and Cq'^ to have non- analytic dependence on m^r. Ideally, all 
dependence should be explicit in the Lagrangian and the coupling Co should only depend on 
short distance scales. However, the fine tuning that results in the large scattering lengths is 
a consequence of a delicate cancellation between long and short distance contributions, and 
in order to put the pole in the physical location, one must induce explicit dependence in 
the perturbative parts of Co [ll6| , |30[| . Eq. ( pTSD will be applied to both S-wave channels. After 
imposing these conditions, there is one free parameter at NLO and two free parameters at 
NNLO. 



III. AMPLITUDES AND PHASE SHIFTS 
A. ^Sq channel 

In this section, we present the NNLO calculation of the ^Sq phase shift. At NLO the 
amplitude involves the diagrams in Fig. |I] calculated in Ref. Q. Graphs contributing to the 
NNLO amplitude include those with one insertion of an order Q operator and those with 
two insertions of either a potential pion or order operator. These graphs are shown in 
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Fig.||. A discussion of the techniques used to evaluate these graphs and exphcit expressions 
for each individual graph are given in Appendix |B[ The NNLO amplitude also receives 
contributions from graphs with radiation pions which are discussed in Appendix |C[ 

By expanding exp{2i6) = 1 + ipMA/{2TT) in powers of Q we obtain expressions for 
the ^5*0 phase shift S = 5^'^^ + S^^^ + S^'^^ (where 5*^"-' ~ Q") in terms of the amplitudes 
A = A-i 

pM Ao 



Ao + Ai (where An ~ Q"), 



1 , / ipM ^ 
In ( 1 + -^A-i 



(2) 



2i 

pM 



271 



Arc 1 + !£M^_, 



^1 



. fpMy 



2n 
2 



Our final result for the amplitude at NNLO is quite simple: 
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Mg\ ml 
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1_ 
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ln(l + 



—m 



An ^87r/2/ Ap 

2 



I 2(f -p2)ImLi2( 



71 



Ap 



2 
2 



2ip 



— ReLi2(^ 

+^ln^fl + ^l-tan-r^hnfl + ^ 



Im Li2 



2ip 
niT, + 2ip 
—rriT, + 2ip 



m 



mt 



(14) 



(15) 



Using Eq. (0) it is easy to verify that the S-matrix is unitary to the order we are working. 
The six linearly independent constants appearing in the amplitude are 7, 0, C2, Cs? C4) Cs- 



7 
C2 

C3 



MCo 
' D2 

{Co? 4/2V47r 



Ci 



Co 



{Co? 



9a Mm^ 
2/2 Air 

27 Mgl 



(Co) 



+ 
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+ 



mt 



(1) 



(0) 



(Co; 



(C, 



(0)n2 



ml STrf^ 



m 

{Co? ' 2/2 UirJ 
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(Co)^ 
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(Co)^ 
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-/^^ + f ) 



M7/M^^2 
47r V87r/2 



ln2-^ 
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^2 



(Co) 



(Co)' 



^2 ZLy2^0 



(Co) 



(Co 



P 47r {Coy 



+ m^ 



P 47r (Co)2 
Ei 2 C2 D2 



rad 



(16) 



C5 



C4 



(Co) 



{C2? ' 

(Co)^ 



Ci — Cs are dimensionless constants. Note that C2 — Cs include factors of m.„ and are not simply 
short distance quantities. After solving the RGB's in Eq. one finds that all quantities 
in square and curly brackets are separately /i independent. Furthermore, the quantities in 
curly brackets vanish at NNLO in the Q expansion due to Eqs. (|Tl|) and (p!^). In Eq. (|T6|) 
the order Q radiation pion contributions appear in given in Eq. ( p37| ) of Appendix ^ 
At order Q, the effect of radiation pions turns out to be indistinguishable from corrections 
coming from contact interactions. 

For the ^5*0 channel, the location of the pole is determined by solving 



0, 



(17) 



This fixes 7 = —7.88 MeV. Note that adding the shape parameter correction to Eq. (|T 
changes the location of the pole by less than 0.01%. The NLO good fit condition in Eq. (|13| 
relates the constants Ci and ^2, 

„2 A/f ^2 /i^ / 27 



ml A-K Slip 



V m^J 



(18) 



leaving one new parameter in the fit at NLO. At NNLO, Cs = once we impose C4 = C|/Co. 
This leaves Cs and ^4, which are related by the NNLO good fit condition 



7 



mi 



2 fMg'j,\2M ml 



^87rj 



ReLi( 



-m^ 



+ 



TT 



12 



(19) 



+ 27^ 

1^) and (19) these conditions basically fix 



47]- ^ 

Since Ci and C4 are multiplied by ^'^ jm?^ in Eqs 
the values of and Cs- We have chosen to fix and (^4 by performing a weighted least 
squares fit to the Nijmegen partial wave analysis pO[. The ranges p = 7 — 80 MeV and 
p = 7 — 200 MeV were used at NLO and NNLO respectively, with low momentum weighted 
more heavily. Using M = 939 MeV, m^ = 137 MeV, = 1-25, and / = 131 MeV the 
parameters for the ^So channel are: 



NLO : 
NNLO : 



Ci = 0.216; 
= 0.0777; 



C2 = 0.0318; 
C2 = 0.0313; 



C3 = 0.1831; (4 = 0.245 



(20) 
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FIG. 3. Fit to the ^5*0 phase shift 5. The sohd hne is the Nijmegen fit |2C| to the data. In a), 
the long dashed, short dashed, and dotted lines are the LO, NLO, and NNLO results respectively. 
In b) we show two other NNLO fits with a different choice of parameters as described in the text. 

The value of these parameters depend on the range of momentum used in the fit, for instance 
using the range p = 7 — 150 MeV at NLO gives Ci = 0.25. From the power counting we 
expect Ci ~ M/(47rA) at NLO and Ci + C4 ~ M/(47rA) at NNLO. For A ~ 300 MeV, 
M/(47rA) ~ 0.25 in reasonable agreement with the fits. 

The ^5*0 phase shift is shown in Fig. §a. The solid line is the result of the Nijmegen 
phase shift analysis [^. The ^Sq phase shift has an expansion in powers of Q, and we plot 
the LO, NLO and NNLO results. The LO phase shift at p = is off by 48%. At NLO, 
the error is 17%. At NNLO, the error in the ^Sq channel is less than 1% at p = m^, and the 
NNLO result gives improved agreement with the data even at p ~ 400 MeV. 

Note that ~ Q is larger than ^2 ~ because from Eqs. ([T8| ) and (|l^), C3/C2 ~ 
m^/(7A7VAr). The parameter ^2 is stable because it is fixed by the NLO good fit condition. 
On the other hand, Ci changes by a factor of 2.8 going from NLO to NNLO. One expects 
the value of coupling constants to change somewhat at each order in the expansion, but a 
factor of three difference is surprising. It is also disturbing that (4 is greater than (i, since, 
on the basis of the RGE and KSW power counting, it is expected that (4 < (i [^. It is 
possible to do a fit and impose the constraints that (i is close to its NLO value and C4 ^ Ci- 
If this is done the error at p ~ is ~ 10%, which is still an improvement relative to the 
NLO calculation and consistent with an expansion parameter of order 1/2. This fit is shown 
as the dotted line in Fig. ^d. 
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The potential diagrams for the ^5*0 phase shift at NNLO were also computed by Rupak 



and Shoresh [16|. To fit (^i and ^4 they essentially demand that the experimental value of the 
effective range is reproduced at both NLO and NNLO. For the observable sin^ 5 at p ~ mjr, 
they find ~ 80%, ~ 65%, and ~ 5% errors at LO, NLO, and NNLO respectively . 

Kaplan and Steele \W\ have proposed that when the perturbative expansion of coupling 
constants is made the sub-dominant couplings should not be treated as new parameters. As 
an example, in their fitting procedure, C2 is given by 

Imposing this condition fixes the value of (^4 so that there is one less free parameter at 
NNLO. Kaplan and Steele motivated this fitting procedure by arguing that adding pions 
should only change long distance physics. Therefore, the number of free parameters in the 
theory with pions should be the same as in the pionless theory. It is worth pointing out 
that Ref. made use of toy models in which the pions were represented as a contribution 
to the potential which is either a delta-shell removed a finite distance from the origin or a 
pure Yukawa. In these models it makes sense to think of the "pion" as purely long-distance 
because the pion effects are cleanly seperated even in the presence of loop corrections. 

In a realistic effective field theory ultraviolet divergences from loops with pions do not 
allow a clean separation of long and short distance scales. As an example consider €2^^'' = 
C2 + C2 + C^^ + . . .. Here C2 first appears in the NLO diagrams in Fig. and introduces 
a short distance effective range-like constant. At NNLO the diagram in Fig. ^ appears and 
has a logarithmic ultraviolet divergence that must be absorbed by C2 ■ This induces a 
ln{fi/K) dependence into the coupling C2 (as is clear from Eq. (pSD ). Since the constant 
K is undetermined it is clear that C2 cannot be determined from lower order couplings. 
Note that if this ln(/i) is instead absorbed into the leading order C2 then this would induce 
additional ln(/i) dependence into the part of the NNLO amplitude that depends on C2. 

In the ^5*0 channel 6*2"^'' does not receive a logarithmic renormaliztion. However, there 
is a new logarithmic divergence that must be absorbed into 02^^ coming from Fig. ^ . 
Therefore, 6*2°'' must be treated as a parameter. It is not possible to renormalize the theory in 
a /i independent way without introducing more parameters than exist in the pionless theory. 
The power law sensitivity to the choice of fi makes the fi independence of observables an 
essential criteria. Since, in general, higher order terms in the expansion of couplings receive 
ultraviolet renormalizations, we prefer to treat all 6*2^'' as free parameters whose size is only 
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^^0 



V2 



V2, 



Fit ^1 
NLO 
NNLO 



2.73 fm 
2.65 fm 
2.63 fm 



-0.48 fm^ 
-3.3 fm^ 
-1.2 fm^ 



3.8 fm^ 
19 fm^ 

2.9 fm^ 



-117 fm^ 
-0.7 fm'^ 



TABLE I. Predictions for terms in the ^Sq effective range expansion. 

restricted by their RGE's. This then imphes that (^4 is a free parameter in both the ^5*0 and 
'^Si channels. However, in the ^5*0 channel at NNLO imposing the relation in Eq. (^) does 
give a yU independent amplitude. In this case the result of the fit is shown by the dashed line 
in Fig. ^3. In general the choice of fit parameters is somewhat arbitrary, and a true test of 
the values can only be made by using them to predict an independent observable. 

Finally, we present NNLO corrections to the higher order terms in the effective range 
expansion 

p COt(5) = + + V2P'^ + f + f + • • • • (22) 



Using the NLO expression for p cot(5), Cohen and Hansen |3T]| obtained predictions for f 2, 



and f4. At NLO, the effective field theory predictions for f2, fs, and disagree with the 
Vi obtained from a fit to the Nijmegen phase shift analysis. The NNLO predictions for the 
shape parameters are shown in Table [HI A| . The prediction for tq is not better at NNLO than 
at NLO, but is still well within the expected errors. The NNLO Vi predictions depend on 
Ci and ^2- We see that the NNLO correction substantially reduces the discrepancy between 
the effective field theory prediction and the fit to the Nijmegen phase shift analysis, but 
the discrepancy is still quite large. This gives some evidence that the EFT expansion is 
converging on the true values of the w,, albeit slowly. Effective field theory predictions for 
the shape parameters have been studied in toy models where one is able to go to very high 
orders in the Q expansion |3^. In the toy models, the effective field theory did eventually 
reproduce the shape parameters, but the observed convergence is rather slow. 
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B. channel 



The S matrix for the ^5*1 and ^Di channels is 2 x 2 and can be parameterized using the 



convention in Ref. [03 



S 



iMp 

271 



^SD jDD 



\ ( e2*^ocos2ei i e"'^''^''^^ si\i2ti^ 



I e 



sin2ei e2*^2cos2ei 
The phase shifts and mixing angle are expanded in powers of Q/K: 



So = + + 5f + . 



6, = + 5t ) +6? + ..., 61 = + 6-« + + 



(23) 



(24) 



The phase shifts and mixing angles start at one higher order in Q than the amplitudes 
because of the factor of p in Eq. (p3D . In the PDS scheme, expressions for 5i^'^\ 6i^\ and 
ei^^ are given in Ref. P]. The prediction for ef"* is given in Ref. [|l9l and is discussed in 

—(2) 

section IIIC, and the prediction for ^2 is given in section HID. 



Expressions for the terms in Eq. (p^) in terms of the scattering amplitude are obtained 
by expanding both sides of Eq. (p3|) in powers of Q. This give^ 



(0) 



2i \ 271 



.(1) _ PM 

"^0 



ASS 
•^0 



47r 1 + ^^^f 



K2) 



pM At 



ss 



47r 1 + ^^^f 

ZTT 



/pMV 



aSS 
•^0 



1 + ^^^f 



+ 



1 + ^^^f 



(25) 



In 5q the terms that depend on Aq^ and Aq^ are purely imaginary and cancel the imag- 
inary part of the term proportional to Af^ as required by unitarity. The order mixing 
amplitude is M 



aSs( 



, +(^ + ^)tan-(^) + ^-2 
^ 37m^N^ ^ ^ 4p^ 



■( 



7m; 



+ 



In 1 + 



The diagrams which contribute to the ^Si amphtude up to NNLO are shown in Figs. |l] and 
and give 



^The branch cut in the logarithm in Eq. is taken to be on the positive real axis. This is 
consistent with 6q{p —>■ 0) = 7c. The sign of our ^Di state is the opposite of Ref. 0, making 
Aq^ in Eq. (^) have the opposite overall sign. 
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The six linearly independent constants appearing in Eq. ( jn\) are: 
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FIG. 4. The '^Si phase shift for NN scattering. The solid line is the Nijmegen multi-energy fit 
|20|, the long dashed line is the LO effective field theory result, the short dashed line is the NLO 
result, and the dotted line is the NNLO result. The dash-dotted line shows the result of including 
the parameter ^5 which is higher order in the power counting. 
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Solving the beta functions in Eq. (P) perturbatively, we find that the quantities in the square 
and curly brackets are separately /i independent, and the quantities in curly brackets vanish 
at NNLO. Q"''^ includes the radiation pion contributions to the amplitude. The expression 
for Q"''^ in the ^Si channel is obtained from Eq. (|C37| ) by interchanging the spin singlet and 
spin triplet labels. 

The LO amphtude ^:^f has a pole at p = ^7 corresponding to the deuteron bound state. 
The deuteron has binding energy B = 2.22 MeV, so 7 = V MB = 45.7 MeV. The remaining 
coefficients, Ci — C4 ci-re fixed using the same procedure as in the ^Sq channel: 



NLO : 
NNLO : 



Ci 
Ci 



0.327; C2 = -0.0936; 
0.432; C2 = -0.0818; 



(29) 



C3 = 0.165; C4 



0.399; 



The ^5*1 phase shift is shown in Fig. ^. The LO phase shift (long dashed curve) has no free 
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parameters, and at p = m^r the error is 60%. The NLO phase shift (short dashed curve) 
has one free parameter {(i), which is fit to the Nijmegen muhi-energy fit (sohd curve). The 
NLO fit to the data is excellent. However, this agreement is clearly fortuitous because the 
NNLO phase shift (dotted line) with two free parameters (Ci, C4) does worse at fitting the 
data than the NLO phase shift. At p = m,r the error is 30%, exceeding expectations based 
on an expansion in 1/2. The error is even greater for larger values of p. The dash-dotted 
line in Fig. § shows the result of including the N'^LO parameter = 0.26. Better agreement 
with the data is found, however, including at this order is in violation of the KSW power 
counting. 

Large NNLO corrections also show up in predictions for the effective range expansion 
parameters. For example, at NLO the effective theory gives an effective range tq = 2.2 fm, 
which is within 20% of the experimental value, tq = 2.73 fm. At NNLO we find tq = 1.3 fm. 
The NNLO correction to tq includes a large negative non-analytic contribution from the 
diagrams with two potential pions. 

The failure of the EFT at NNLO in the ^5*1 channel is due to large corrections from the 
two pion exchange graphs in Figs. ^,k,m. The term which dominates the NNLO amplitude 
for large p is 

^f''-i^-?'^s©)^'--'(;^)- '») 

For p ^ m^r this term grows linearly with p, an effect which can be clearly seen in Fig. ^ (the 
growth in Fig. ^ is quadratic due to the extra p in Eq. (P5|)). The contribution in Eq. ( |3D[ ) 
is large because of the coefficient of 6 which is much greater than the expansion parameter. 
For p ^ the size of this contribution relative to the LO amplitude is Snp"^ / Aj^^. The 
fact that this correction survives in the chiral limit indicates that it comes from the short 
distance part of potential pion exchange. Large non-analytic corrections are also found in 
some of the other spin triplet channels at this order. 



C. ^Si - channel 

The ^Si — ^Di mixing amplitude at NNLO was presented in Ref. The result 

is briefly summarized here for the sake of completeness. The prediction is shown in 
Fig. 1^. For ei the LO (order Q^) prediction vanishes and the NLO prediction is pa- 
rameter free. At NNLO there is one free parameter C'^^\ij) which is fit to the data: 
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FIG. 5. The "^Si — mixing angle for NN scattering. At LO this phase shift is zero. The 



dotted hne is the NLO result ||9[ and the dash-dotted line is the NNLO result [19|. The solid line 
is from Nijmegen's multi-energy partial wave analysis pO|. 



a 



(5D), 



— 4.6fm^. This value is consistent with the power counting estimate which 



gives |C'2^^'*('^7r)| ~ Air / (Mm'^A) ^ 3.6fm^ for A = 300 MeV. The mixing angle agrees with 
expected errors for p ~ 171,^, but for larger values of momentum there is serious disagreement 
between theory and experiment. For p ~ 171.,^ this disagreement is comparable to the uncer- 
tainty of a calculation of ei within the Weinberg approach ||^. A more recent analysis ^ 
gives a more accurate prediction for ei, but an analysis of the uncertainty due to the cutoff 
dependence is not presented. 

At the order we are working A^^ = Aq^ + ^ and 
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The behavior of this mixing angle for p ^ rrij^ can be examined by taking the m.,^ 
of the mixing amplitude: 
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where the ln(/i) dependence is cancelled by C'^^\ In this channel the term proportional 
to irp is suppressed by an additional factor of 7, and the dominant terms in the NNLO 
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FIG. 6. The order Q° dia gram (a) and order Q diagrams (b and c) that contribute to the P 
and D wave channels. Only the ^Di channel gets a contribution from diagram c). 

calculation for p ^ m.,^ are analytic, growing as p^. The fit to the low energy data in 
Ref. |]19[ did not give a value of C2^^^ that cancelled this growth as can be seen clearly in 



Fig. |. 

An interesting way to test the EFT for nucleons is to compare the value of Cf^\ij) 
extracted from our NNLO calculation of ei to the C'f^\ji) extracted from the NNLO 
calculation of the deuteron quadrupole moment [|l^]. To make the comparison meaningful 
the same renormalization scheme must be used (and the same finite constants must be 
subtracted along with the /e pole). Ref. [0] does not explicitly give C^^^ counterterms 
so it is was not possible to compare our value for cf^\mT,) with the value extracted there. 



D. ^Di channel 



In the KSW expansion, there is no order 1/Q contribution to the ^Di — '^Di amplitude. 
Using Eq. (E3]) we can derive expressions for the ^Di phase shift up to order Q^: 
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The last two terms in ^2^^ are purely imaginary and cancel the imaginary part of A^^ as 
required by unitarity. The NLO contribution comes entirely from one pion exchange 
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(34) 



Four nucleon operators which mediate transitions between two '^Di-wave states must have at 
least 4 derivatives. Graphs with these operators do not contribute until order in the KSW 
expansion. (The leading operator which mediates '^Di wave transitions is renormalized by 
graphs with two insertions of the O^^^ operator. An insertion of O^^"* is order Q, therefore 
these graphs are order .) At NNLO, the ^Di — ^Di amplitude gets contributions from the 
graphs in Fig. |^. 
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10 20 30 40 50 100 , , 200 300 

p(MeV) P(MeV) 

FIG. 7. ^Di phase shift for NN scattering. The sohd line is from Nijmegen's multi-energy 
partial wave analysis ||2^. At LO this phase shift is zero. The long dashed line is the NLO result 
and the short dashed line is the NNLO result. In the first plot the NNLO prediction without 
including the graph Fig. |^ is displayed as a dotted line. The only parameter entering at this order 
is Cg which is fixed by the location of the pole in the ^5*1 amplitude. 



(■3 e \ 

The only short distance operator which contributes to this amplitude at NNLO is i 
whose coefficient is completely determined by the location of the pole in the spin-triplet 
channel. Therefore, no free parameters appear in the calculation of this amplitude. The 
NNLO amplitude is: 
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Values for the individual graphs are given in Eqs. (P27|) and ( P29|) . 

The NLO and NNLO predictions for 82 are plotted in Fig. 0, along with the result of the 
Nijmegen partial wave analysis. The NLO result gives satisfactory agreement with data up 
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to 300 MeV. The NNLO calculation is less accurate than the NLO calculation especially for 
p > 50 MeV. The error in the NNLO calculation is always greater than the NLO calculation, 
so for this observable there is no sign of convergence of the KSW expansion at any value of 
p. At NNLO the prediction for the ^Di phase shift suffers from the same problem as the ^5*1 
phase shift, namely a large term in the amplitude that grows linearly with p for p ^ rriT,-. 
Taking m,r ^ we find 

^ ADD _ -1 I ^ f^P , 37rpx 

4vr'^ 2A^^ + A^^ U + 2(7 + ^p) 14 ) ' ^^^^ 

The last term in this equation dominates the phase shift at large momenta. 

Note that for low momentum, the inclusion of graph c) in Fig. |^ improves the agreement 
over a theory which contains only perturbative pion exchange. This can be seen in Fig. |^ 
where the small dashed line (NNLO with c)) lies closer to the Nijmegen phase shift (solid) 
than the dotted line (NNLO without c)). 

In this section we have presented calculations of the phase shifts and mixing angles in 
the ^5*0, ^5*1, and ^Di channels at NNLO. We found that the ^5*0 phase shift agrees well with 
data up to p ~ 400 MeV. However, in the spin triplet channels the effective field theory 
expansion does not seem to converge. The ^5*1 — ^Di mixing angle ei agrees with data 
to within errors for p ^ m^r. This is not true for the ^Si and ^Di phase shifts. In these 
channels, two pion exchange graphs give corrections which worsen the agreement with data. 
This suggests that the perturbative treatment of pions is inadequate in spin triplet channels. 



IV. P AND D WAVE CHANNELS 

In this section we will examine the ^Pi, ^Po,i,2, ^-02, and ^-D2,3 channels in an effort 
to get a better understanding of perturbative pions. In these channels there is no order 
contribution, the contribution consists solely of single pion exchange (Fig. ^ja), and 
the order Q contribution comes from the potential box diagram in Fig. ^o. Four nucleon 
operators only contribute at higher orders in Q. Since the coefficients of these operators are 
not enhanced by the renormalization group flow near the fixed point they have a scaling 
determined by dimensional analysis. In the P waves contact interactions first appear at 
order Q^, while in the and "^-02,3 they first appear at order Q^. 

In Ref. []rT| the phase shifts with L > 2 were calculated using perturbative pion exchange. 
In this calculation, the one loop potential box, soft diagrams, and a subset of order 
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FIG. 8. Order soft pion loop graphs for the P and D wave channels. 

corrections were included simultaneously. The potential box in Fig. |^ is order Q, while the 
soft diagrams in Fig. |^ are order Q^. At order there are also relativistic corrections and 
radiation pion contributions. The latter can be absorbed by using the physical value of qa 



in the one-pion exchange diagram. However, Ref. fTj] did not include the double potential 
box 

lAVi^^%^, (37) 



I I I 
I I I 

J I L 



V2/2/ 

which is also order Q^. Since a complete order amphtude is not yet available, no diagrams 
of order or higher will be included in our analysis. 

The order Q, phase shifts are given in terms of the amplitude by: 

5« = ^4^ ^(^) = ^Ref41. (38) 

Att ' 47r L M ^ ^ 

Projecting Fig. |^a onto the various P and D waves using the projection technique discussed 

in Appendix ^ gives the results in Eq. ( |A9| ) which agree with Ref. |Tl[]. In these channels. 



the box graph in Fig. |]b can be evaluated analytically using the techniques discussed in 
Appendix B. We have instead chosen to calculate the partial wave amplitudes by using 



the expression for the box graph given in Ref. [|TT|, and doing the final angular integration 
numerically. 

Results for the P and D wave amplitudes are given in Figs. ^ and |10| respectively. The 
potential box diagram gives a very small contribution in the singlet channels, in contrast 
to the triplet channels. For momenta p ^ 110 MeV the NNLO calculation gives reasonable 
agreement in the ^Pi, ^Pi, ^D2, and '^-02,3 channels, but not in the '^Po,2 channels. For larger 
momentum, p ~ 300 MeV, the error in the ^Pi, ^Pq^i, and channels is very large. This 
is less of a concern because the KSW power counting is not designed to work for p ^ m^r, 
but it does indicate a need to modify the KSW power counting for momentum greater than 
the pion mass. 

To get a better idea of what is happening at large momenta it is useful to look at the 
rriT^ —>■ limit of the amplitudes: 
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FIG. 9. P-wave phase shifts for NN scattering. The triangles are values from Nijmegen's single 
energy analysis and have errors that are invisible on the scale shown. At LO these phase shifts are 
zero. The long dashed line is the NLO result and the short dashed line is the NNLO result. There 
are no free parameters at this order. 
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FIG. 10. D-wave phase shifts for NN scattering. The triangles are values from Nijmegen's 
single energy analysis and have errors that are invisible on the scale shown. At LO these phase 
shifts are zero. The long dashed line is the NLO result and the short dashed line is the NNLO 
result. There are no free parameters at this order. 
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TABLE II. Comparison of the NNLO part of the phase shift to its niTj — > hmit at 
p = 300 MeV. 

47r^^ - 7 Ann + A^^ + 49 ^ 

For the spin singlet channels there are no corrections which grow with p, while the spin triplet 
channels have non-analytic corrections proportional to np. This short distance behavior is 
similar to what is seen in Section III. At p = 300 MeV these particular non-analytic terms 
dominate all other NNLO corrections as can be seen from Table 0. In the ^P2, "^-02,3 channels 
these corrections improve the agreement with data, while in the ^Po,i channels they do not. 

At lower momenta p ~ 50 MeV the effective theory does a better job of reproducing the 
phase shifts. Therefore, it seems possible that in these channels predictions for terms in 
the effective range expansions, p'^^~^^ cot 6 = — 1/a + + . . ., might work fairly well. 

Equivalently one can match onto the theory with pions integrated out to make predictions 
for the coefficients of four nucleon operators in the P and D waves. Such an investigation is 
beyond the scope of this paper. 



V. DISCUSSION 

In this section we summarize our results for the S, P, and D wave phase shifts. We also 
discuss in greater detail the nature of the perturbative expansion in the spin singlet and 
triplet channels. 



Errors in each channel at p = 50 MeV and p = are given in Table |T|. For an 
expansion parameter of 1/2, we expect roughly 50% error at NLO and 25% error at 

NNLO [Q). (For the two S-wave phase shifts, which start at one lower order in the expansion, 
the expected error at NLO and NNLO is 25% and 12.5%, respectively.) At p = m^r, errors 
are significantly larger than expected in the ^5*1, 'P02, and ^Di channels. In the case of the 



3 



D3 channel the percent error is exaggerated due to the smallness of the phase shift, so in 
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p= 137 MeV 
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^Pi 
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NLO 
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33% 
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TABLE in. Percent errors in the phase shifts for p = 50 MeV and p = 137 MeV at NLO and 
NNLO. (* Since the ^1)3 phase shift is cfose to zero the percent errors are not very meaningful. 
The absolute errors for this phase shift have the expected size.) 



this channel the percent error is probably not a figure of merit for examining the quality of 
the expansion. However, the LO correction in this channel has the wrong sign so there is 
no sign of convergence of the perturbative expansion. At p ~ m^r, the performance of the 
effective theory is erratic, working some but not all of the time. The overall agreement with 
data at 50 MeV is better, but there are still channels (^Po,2; '^-Di) in which the agreement 
with data is worse than one expects. 

The perturbative expansions in the spin triplet and singlet channels are qualitatively 
different. All triplet channels have non-analytic corrections that grow with p, while the 
singlet channels do not. This can be understood as follows. First consider the spin singlet 
channel. In this channel, the potential due to one pion exchange is the sum of a delta 
function and a Yukawa potential, 

g'2 _|_ ^2 (^^2 _|_ ^2 -) 

The effect of the delta function part of one pion exchange is indistinguishable from the 
Cq ° operator and therefore only contributes to S-wave scattering. A well known theorem 
from quantum mechanics shows that at large energy the Born approximation becomes more 
accurate for a Yukawa potential. From the point of view of the field theory, this means that 
the ladder graphs shown in Fig. ^ with Yukawa exchange at the rungs are suppressed by 
powers of the momentum. Using dimensional analysis we see that adding a Yukawa rung 
gives a factor of 
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FIG. 11. Potential pion ladder diagrams. 
g\M r (T-k ml ml 



(41) 



2/2 y (27r)"(A;2 + 2fc-p) (A;2 + m2)... p Ann 

for p ^ m^r. The one loop pion box diagram in the spin singlet channel gives a contribution 
which can be eliminated by a shift in Cq " and other terms that are suppressed by powers 
of m^/p (see for e.g., Eq. ( P18|) ). Once the short distance effects of pions are absorbed into 

/lex 

C'o the remaining piece of two pion exchange is never larger than the estimate given in 
Eq. ([^) and gets smaller as p increases. This is good for the convergence of the perturbative 
expansion because it means higher order potential pion corrections in singlet channels will 
be well behaved. In fact the tree level pion exchange graph gives almost the same prediction 
for the ^Pi and phase shifts as the leading order prediction in the Weinberg expansion 
that sums potential pions to all orders. Thus, the evidence for the behavior of the spin 
singlet channels is independent of how the parameters in the ^5*0 channel are fit to the data. 

The S, P and D wave phase shifts are calculated to NNLO within the Weinberg expansion 
in Refs. and These studies are complementary since Ref. gives the uncertainty 
in their NNLO predictions by varying the cutoff from 0.5 GeV to 1.0 GeV, while Ref. 
explicitly displays their LO, NLO and NNLO results (which are respectively order p"^, 
and p^ in the potential). In comparing our results with those of Refs. [^,0 it must be 
noted that these calculations include many effects which do not appear until higher order 
in the KSW expansion. For example, P-wave contact interactions are included at NLO so 
the predictions in these channels have a free parameter which is fit to the data. Soft pion 
effects^ also enter at this order. These effects enter at order in the KSW power counting 
(N^LO). In the Weinberg expansion the singlet and some triplet phase shifts cannot be fit 
until these interactions are included [0,0. 

At p = 306 MeV the LO result in Ref. is 6{^Pi) = -12° and 6{^D2) = 2.2° which is 
very close to tree level pion exchange which gives S{^Pi) = —12° and 6{^D2) = 2.1°. Thus, as 



^Soft pion diagrams with nucleons and A's were calculated in Refs. [pT] , ^ 
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TABLE IV. Predictions for the phase shifts when diagrams with insertions of Cq and potential 
pions are summed to aU orders. The results shown are the leading order predictions from Ref. [^. 

expected, the discrepancy between theory and experiment seen in the ^Pi and channels 
in Figs. P and jl^ is not removed by summing potential pion diagrams. The LO predictions 
in the Weinberg expansion are shown in Table In the ^ 5*0 channel the result in Table 
is only slightly better than the LO result in Fig. We conclude that there is little to be 
gained by summing potential pions in spin-singlet channels. 

In the spin triplet channel, the potential from one pion exchange is much more singular 
and has terms that go like for small r (where r is the nucleon separation). In fact, 
without introducing an ultraviolet cutoff, it is not possible to solve the Schrodinger equation 
for such a singular potential. In field theory this means that higher potential pion ladder 
graphs have ultraviolet divergences of the form p'^"^/e, which must be cancelled by a four 
nucleon operator with 2m derivatives. (Examples of two loop graphs with p'^/e divergences 
were computed in Ref. [^.) In the spin triplet channel perturbative pions give corrections 
which go like {p/Ann)^, where k is the number of loops. Loop graphs with pions in the spin 
triplet channel can therefore have finite corrections which grow with p and are non-analytic 
in p^. These short distance pion corrections cannot be compensated by any short distance 
operator. For the ^Si, ^Pq, and ^Di channels the non- analytic corrections are large and ruin 
the agreement with the data. In the ^P2 and ^D^ channels the quality of the expansion is 
poor because the non-analytic correction makes the NLO and NNLO corrections comparable 
in size. 

In the ^Si channel calculations within the Weinberg approach [^,0 have small cutoff 
dependence and agree much better with the Nijmegen partial wave analysis than Fig. ^. 
The ^Pi and ^Di channels are also in good agreement. In these channels the summation 
of potential pions improves the agreement with data. However, in other P and D wave 
channels the summation of potential pions is not as helpful. At LO Ref. [0 finds large 
disagreement in the ^Pq and ^P2 channels as can be seen from Table |rV|. These predictions 
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are similar to what is given by tree level pion exchange. At NNLO the ^Pq phase shift is 
in reasonable agreement with the data with small cutoff dependence 0. In the channel 
there is larger cutoff dependence and no sign of convergence of the perturbative expansion 
at p ~ 300 MeV 0. In the ^Po,2 channels agreement with data is only achieved at the order 
that a free parameter appears. Soft pion graphs and four nucleon operators appear to be 
more important than summing potential pions. In the '^-02,3 channels at p ~ 300 MeV our 
NNLO prediction is of similar quality to the NNLO prediction in Ref. 0, so the summation 
of potential pions does not seem to be necessary. 

The large NNLO corrections in the ^Si, ^Di, and ^Po,2 channels cast considerable doubt 
on the effectiveness of the KSW power counting for pions. The ~ 10% accuracy of NLO 



results [|T2| remains somewhat mysterious. For momenta p < m^r the pion can be inte- 
grated out. This low energy theory has been shown to be effective in calculations at NNLO 
|2|,|3|^. 



Large perturbative corrections from two pion exchange suggest that a nonperturbative 
treatment of pions is necessary for nuclear two body problems in some spin triplet channels. 
This is achieved in Weinberg's power counting because the potential pion exchange di- 
agrams are summed at leading order. However, the graphs which are resumed by solving 
the Schrodinger equation have logarithmic divergences of the form m^™ ln(A) or ln(A) 
(this has been shown explicitly for the case m = 1 in Refs. ||^ and |^). The short distance 
counterterms necessary to cancel the ln(A) dependence of these graphs are not included until 
higher order (for m = 1 this would be C2P^ and D2Tn^). The residual cutoff dependence is 
of the same size as higher corrections. However, it is not a priori clear why the contribution 
of the graphs included in the summation is larger than the omitted counterterms (see for 



instance Ref. ||37|1). 

We have seen that in spin triplet channels there are perturbative corrections which sur- 
vive in the m.,^ —>■ limit, and are enhanced by large numerical factors (the vr's in Eq. (|39D). 
Empirically these terms tend to dominate the NNLO correction. Furthermore, these cor- 
rections are non-analytic in p^. Since no unknown counterterm can contribute to their 
coefficient they can be calculated unambiguously. It would be interesting to see if there 
are similar large calculable corrections at higher orders. If so, then the advantage of the 
Weinberg approach relative to the method used here is that it sums these important con- 
tributions along with smaller scheme dependent corrections. In three body problems a 
power counting similar to KSW gives accurate results at very low energies. In these com- 
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putations, the perturbative treatment of pions and higher derivative operators is crucial 



because it renders the calculations more tractable (see Ref. [39|) than conventional potential 
model approaches. For this reason, an approach to two body forces which sums genuinely 
large calculable corrections from pion exchange analytically or semi-analytically would be 
worth pursuing. 

In this paper we extended calculations of the ^5*0, '^S'l, and ^Di phase shifts to NNLO in 
the KSW power counting, including a complete calculation of radiation pion contributions. 
At this order the predictions for the P wave and remaining D wave channels were also 
examined. In spin singlet channels a perturbative treatment of potential pions is justified. 
The large disagreement for the ^5*1 phase shift provides an unambiguous indication that the 
KSW expansion for pions needs to be modified. This is supported by the failure in the ^Pq 
channel and the lack of convergence at p ~ in the '^-Di, ^P2i and '^D^ channels. 
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for their hospitality. T.M and I.W.S. were supported in part by the Department of Energy 
under grant numbers DE-FG03-92-ER 40701 and DOE-FG03-97ER40546. 
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APPENDIX A: PARTIAL WAVE PROJECTION TECHNIQUE 



In this Appendix we discuss a method for obtaining the contribution of a Feynman 
diagram to a particular partial wave amplitude. We use a trace formalism which allows us 
to project out the partial wave amplitude before doing loop integrations. This approach has 
the advantage of being well adapted to situations in which spin (and isospin) traces should 
be performed in n = 3 — 2e dimensions. 

Consider the process N{k/2 - p)N{k/2 + p) N{k/2 + p' )N{k/2 -p'). We begin by 
defining two nucleon states |28[1 

1 



NN{s;k,p) 



p 



N^{k/2 + p)P^'^N{k/2-p) |0), 



(Al) 



where s = "^^^^Lj, p = \p\, and the matrix pi^^'^^^j) projects onto the desired partial wave. 
The normalization of the states in Eq (|A1|) is chosen so that averaging over polarizations 

(A2) 

with the projection matrices satisfying 



(NN{s';k',p') \NN{s;k,p)) = 6%k' -k)6{p' -p)6''' 

pol. avg 



(A3) 



pol. avg 

Here Tr denotes a trace over spin and isospin. Evaluating the traces in n dimensions gives 
the following normalization to the projection matrices for the S, P, and D waves: 

pi's,) ^ (^^2 ^ ■ e) {ir2) 



(i5o) _ ii(^2) iir2T ■ e,) 
2V2 
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2V2 



(A4) 
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where e, e^^ , and 6*-^'^ are the J = 1, 2, 3 polarization tensors and p = The ^5*1, ^Pi and 

^-Di,2,3 waves are isosinglets, while the ^Sq, ^-Po,i,2 and states are isovectors labelled by 
ej. Averaging over the polarization states in n dimensions gives 



1 



n 



n{n - l){n + 4) In + 2 



{n + 2){n-l) 
-2 



n 



(A5) 



+ 



(^gd^jm^kq ^ ^il^jq^krr,^ ^ (^i ^ j ^ k ^ i) + {i 



J 



To evaluate the matrix element of an operator O, we write O = N*N^Oab-cdNcNd, so the 
scattering amplitude is 

dcosO 



iA= J2 {NN{s']k',p)\o\NN{s]k,p)) = A j ^ 



Pab ^ab : cd Pcd ' 



(A6) 



pol. avg 

where p' ■ p = p^ cos 9 and the indices {a,b,c,d) are for both spin and isospin. 
Examples of the use of Eq. ( [A6| ) are: 



CSo) 



'So X 'So = (- / rfcos^^) 4 (-zCo)Tr 



pCSo)peSo) tl rj.[pCSo)p{^So) t 



-iCo , 



(AT) 



where ' is given in Eq. (H) and we have averaged over the isospin polarizations, and 



_j_ I ^ ^ Tr 



p(^5o) ^. _ ^ ^kpi'So) t ^T. _ ^) 



{p' — p}"^ + ml 



(A8) 



where in evaluating this trace it is useful to recall that (io"2) cr^ (^o"2) = <?. The factors of 4 
and 2 in Eqs. ( |A^ ) and ( [A8| ) are symmetry factors for the graphs. Projecting the tree level 
one pion exchange diagram onto the various P and D waves gives the order amplitude 
in these channels: 
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These expressions agree with Ref. JTT 



APPENDIX B: EVALUATION OF ORDER Q LOOP DIAGRAMS 

In this Appendix exphcit expressions are given for the individual graphs in Fig. ^ in the 
^5*0 and ^5*1 channels and the graphs in Fig. |^ for the ^Di channel. Details on the evaluation 
of the three non-trivial two pion exchange diagrams (Fig. |]i,k,m) are also presented. 

Our calculation is performed using the Power Divergence Subtraction (PDS) renor- 
malization scheme in d = n + 1 dimensions. A factor of (;u/2)^^" is included with each 
loop and we work in the center of momentum frame, N{p)N{—p ) N{p' )N{—p' ). A 
detailed description of the method used to implement the PDS scheme can be found in 
Ref. Our results are slightly different than Ref. because all spin and isospin traces 



are performed in n dimensions rather than 3 dimensions. For a four-nucleon operator with 
coupling C, there are subtractions for ultraviolet divergences in n = 3, 5"^C, and we define 
the renormalized coupling C(/i) by: 

oo 

(^bare _ ^finite _ ^uv^ ^ ^finite _ _ ^ ^"^^'(/i) . (Bl) 

m=l 

Here 6"^C{fi) is the finite m-loop PDS counterterm, which is defined by canceling overall 
poles in n = 2 (linear divergences) and then continuing back to n = 3. This procedure 
correctly accounts for the unusual scaling of the four nucleon operators due to the presence 
of the non-trivial fixed point. C(/i) may also cancel ln(/i) dependence in the amplitude. The 
beta functions in Eqs. (H) and (|) are computed using 

P = f^^Cif^) = E I- . (B2) 

Renormalized PDS diagrams are defined by adding graphs with counterterm vertices to the 
original diagram. 
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FIG. 12. Reduction in the topology of non-relativistic loop graphs from performing the energy 
contour integrals and picking poles from the marked nucleon propagators. 



1. Basic Strategy for evaluating non-relativistic loop integrals 

Our basic strategy for evaluating massive multiloop potential diagrams analytically con- 
sists of the following three steps: 

1. Evaluate the spin and isospin traces, then do the energy integrals using contour in- 
tegration. This leaves integrations over loop three-momenta which will be evaluated 
using dimensional regularization in = 3 — 2e dimensions. When nucleon poles are 
taken in doing the contour integrals in the n + 1 dimensional non-relativistic theory, 
the remaining loop integrals have the same form as n dimensional loop integrals in 
a Euclidean relativistic theory. The corresponding diagram in the n dimensional Eu- 
clidean theory can be found simply by shrinking to a point the nucleon propagators 
whose pole is taken. This gives a graph with a "reduced topology" . Two examples of 



this are given in Fig. |T2[ In the first example the energy integrals are performed using 
the poles in the marked nucleon lines and the two loop graph becomes a two-point 
function. Only one momentum is relevant to the evaluation of this diagram because in 
the original graph the loops only depend on the relative momentum between the two 
outgoing lines. In the second example choosing nucleon poles as indicated the three 
loop graph becomes a vacuum bubble. In the original diagram the loops only see the 
total incoming energy. This energy will appear in mass terms in the reduced diagram. 

2. Eliminate factors of momenta in the numerator. We begin by canceling terms in the 
numerator against terms in the denominator (partial fractioning). Numerators which 
can not be reduced by partial fractioning are labeled irreducible. These numerators 



are dealt with using the integration by parts technique 40], using the tensor decompo- 



sition technique ||41[], and/or by using relations due to Tarasov [^. Tarasov's method 



is to derive relations between integrals in n dimensions with irreducible numerators 
and integrals in n + 2,n + 4, . . . dimensions with trivial numerators. These integrals 
are then reduced to n dimensional integrals with trivial numerators. (This method 
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was automated for two loop graphs in Ref. [H3| using a Mathematica program called 



Tarcer). A review of these techniques is given in Ref. [H 



3. Evaluate the remaining scalar integrals. This can be done directly using Feynman 
parameters, however it is often more useful to switch to position space using 

^.-^ = / (f^R e-''^-^ GiR, m) , (B3) 

where G is the position space Green's function. An m-loop momentum space integral 
with k propagators becomes a (A; — m)-loop integral in position space. In n dimensions 
the Green's function is 

G[R^m) = [J^^J^ = '('^Y''-'k,.^,MR) , (B4) 

where i^' is a modified Bessel function. For odd n the Bessel function becomes an 
exponential; for n = 3, G{R,m) = e~™^/ (47ri?) . If the reduced topology is that of a 
zero or two point function there are no non-trivial angular integrations. Since these 
are exponential integrals the finite part of graphs are easy to evaluate. To evaluate 
ultraviolet divergent integrals we follow Ref. and split the d"^R spatial integration 



region into two parts /q°° dR = Jq dR+ dR. Ultraviolet divergences occur for i? ^ 
so the dR integral can be done with n = 3, discarding terms that vanish as L — >■ 0. 
For the Jq dR integral we expand the Bessel functions about R = using 

and then do the integration. Ultraviolet divergences are expressed as 1/e poles just as 
if the integration had been carried out in momentum space. When the integrals from 
to L and from L to cxd are added the L dependent terms cancel. For n = 3, the 
scalar two-point and vacuum diagrams with arbitrary masses have been evaluated to 
two and three loops respectively in Ref. 



2. The order Q potential diagrams 

At order Q the potential diagrams that contribute to S-wave NN scattering are shown 
in Fig. |[ The evaluation of the graphs in Fig. ^-h,j,l,n is the same in the ^Sq and ^5*1 
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channels, while Fig. ^,k,m differ. In the ^5*0 channel the order Q diagrams have also been 
evaluated in Ref. [|1^], however our results are slightly different since all traces are performed 
in n dimensions. 

The graphs in Fig. ^-f are simple to evaluate: 

2C2-D2\ 2 2 , n-\A 1 3 ^2-^2 2 2 



h)+e) = -^ [A^rY - j + 2 ^ [A^if ^ p'mi , (B6) 



XCY {Co)y - ' ' (Co) 

The diagrams in Fig. Hg,h are also straightforward. Renormalized diagrams are calculated 
by adding diagrams with the appropriate PDS counterterms. The two basic renormalized 
diagrams needed to evaluate the diagrams in Fig. ||g,h are: 




2% {C^p' + D^mi, 



^"^^ln(l-^) 



2p ^ m. 

><!>< - ^ (^-^ - f f >" (^) 



+ r'^CW,, (B7) 



2, 2ip^-\ .MA-i,. ,Mgl^ . 



where the sum of the C2 and D2 operators is represented by a diamond. The diagram in 
Eq. ( p8l) is ultraviolet divergent and in defining the renormalized graph we have introduced 
two counterterms to cancel the 1/e poles {n = 3 — 2e): 



2 47r 87r/2 

The PDS renormalization scheme is being used, so there are also finite subtractions that 



correspond to poles in three dimensions. The graph in Eq. (B7) does not require a PDS 
counterterm because we are evaluating spin and isospin traces in n dimensions, and the 
isospin traces gives a factor of n — 2 which cancels the l/(n — 2) pole in the loop integration. 
The graph in Eq. (|B^ ) has C2 and D2 PDS counterterms which produce the factor of /i^, 
while the factor of fi is from the first graph with a 6^C2 counterterm in place of the diamond. 
The remaining diagrams in Fig. |^,h follow by dressing the results in Eqs. ( |B7[ ) and ( [B8| ) 
with Co bubbles and adding the appropriate counterterm diagrams. The final result for the 
^5*0 and ^5*1 channel is 
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Next consider the graphs in Fig. § with two potential pions. Diagrams j), 1) and n) can 
be obtained using the expressions for the NLO one pion exchange diagrams: 



2 r 



—I 



1] 



ip 



M[A^,f Mg\ 



^m^.^_2_ip. 



2p 



p -/i 



2 



47r 87r/2 

giving the following expressions valid for the ^5*0 and ^Si channels: 
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The last diagram required a new ultraviolet counterterm 

1 , o , A/T, ^ , ^2 
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while the other poles in the graphs in Eq. ( |B13| ) are cancelled by diagrams with the D2 
counterterm defined in renormalizing the graph in Eq. ( |B12| ). 

To evaluate the diagrams in Fig. ^,k,m we follow the three steps discussed in section 



In the ^Sq channel step 2 may be accomplished by canceling terms in the numerator against 
those in the denominator. For example, after doing the contour integrals the integrand of 
the one-loop box diagram is 

dJ'k k'^{k-qY 



{2nY[k^ + 2k-p][k^ + ml][{k-qY + ml] 
(Fk 1 



(B15) 



ml 



{2nY [k^ + 2k-p] 



ml 



[k'^ + m'l] [{k-qY + ml] [k'^ + ml][{k-qY + m'l]) 

where q = p'—p. The integral over k can be evaluated using Feynman parameters. The term 
with three propagators requires the most effort and gives an answer involving di-logarithms, 
Li2. Integrating over cos6 = p ■ p'/p"^ to project out the ^Sq partial wave gives 
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Manipulations similar to those in Eq. ( B15 ) allow us to eliminate the numerators in Fig. ^ 
and Fig. |^m. For these diagrams all the remaining scalar integrals were evaluated by Rajante, 
in Ref. A D2 counterterm is introduced to cancel an ml/e divergence in Fig. ^|m, 
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The final result for Fig. |^i,k,m in the ^Sq channel is then: 
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(B18) 

-m,r + 22p^ 12 J J 

Only the three loop graph requires a PDS counterterm because the isospin trace with two 
pions gives a factor of [n — 2)' while each loop gives at most a l/(n — 2) pole. Our analytic 
expression for the box diagram agrees numerically with the result in Ref. |Tl[ . 

The evaluation of Fig.^,k,m in the ^5*1 channel is more difficult because of the more 
complicated numerators. For the box graph we can again perform step 2 of the previous 
section by partial fractioning, 

4[(A?-g)-A?]2 + (n-4)A^2(fc-g)2 [ d^k 1 
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2ml " 2^' 



+ 



3m^ + Amlq^ + q 



[(fc — g)2 + m2] [k"^ + ml] [k'^ + ml][{k — qY + ml]} 

Since this graph is finite we have set = 3. The terms with three propagators require 
Eq. ( [B16|) and the following two integrals 



41 



dcos9{l-cos9) J 
1 Uml + 2p2^ 



d^k 1 

(27r)3 [k2 + 2k ■ p][k^ + ml][{k - qf + ml] 



tan 



ml + 2p2 



2p' 



7M^ + 



mi 



4p' 



2p^ — ip 



+ Im Li2(- 



-2p^ + ipmj, 



mi 



) 



/ wi'n , / 4p \ z /ml ml\ , 9 / 4» \ 



J^^dcos9{l-cos9y J 



m 
d^k 



4 V pi 



4p2 



X tan" 



3m^ /15m^ 



m-j^p 
ml + 2p^ 



(27r)3 [A; 2 + 2A; ■ p ] [/t 2 + m2] [{k - g )2 + m2 

3 

2^ 4 

4 



9ml \ . /_ p2 



16p^ 6Ap 



V 2p2A ^ p2 ^ 4„4 



N 9 / ml mZ\2 I 



4p4 



' 2p — ip m^ 
m2 + 4p2 



+ Im Li2 ( 



4p4 

— 2p2 + z p m^ 



^ 3m2 ^ 3ml 



mi \ 2 . 



4p4 



V 2p2 / \ n2 / V 7772 



4p" 
rnt 



+ + ln2(l + -^ 

4 V „2 4„4 ) V rn2 



4 V p2 4p4 

Using these results we find that the renormahzed box graph in the '^S\ channel is 
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The /i dependence comes from adding a S^Cq counterterm at tree level to cancel a l/(n — 2) 
pole. For /i = 0, Eq. ( |B21| ) agrees numerically with the result in Ref. 
For the ^Si channel, the two loop graph in Fig. ^ requires evaluating 
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(271)^^ y (27r)" [^-i + 2l-p'][k^ + 2k-p'][k^ + ml][{I-ky + ml] ' 
We begin by eliminating the loop momenta from the numerator. This may be done using the 
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computer program^ in Ref . |^ , that implements a set of reduction formulae due to Tarasov 
J2|. The remaining scalar integrals can then be found in Ref. |56]. We have checked by 



hand that this program gives the same final result as using tensor decomposition along with 
integration by parts and partial fractioning. The following counterterms are needed to cancel 
l/(n — 3) poles: 
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We find that in the ^Si channel the PDS renormalized diagram is 
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The term proportional to is from a 6^Co{fi) counterterm, while the term proportional to 
/i is from a one- loop nucleon bubble with 6^Co{fi) and Co{fi) vertices. 

Now we turn to the three loop diagram in Fig. |^m in the ^5*1 channel. After performing 
the traces and energy integration we are left with the integral: 

4[ (g - f) ■ (g - k) f + {n- 4) (g - if (g - 
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To eliminate the first term in this numerator we have implemented by hand the procedure 



given in Ref. . The remaining 3 — 2e and 5 — 2e dimensional scalar integrals are evaluated 
in position space as described in Step 3 of the previous section. The non-analytic ultravi- 
olet divergences in the result {p^/e, m^p/e) are cancelled by inserting the counterterms in 



^For this loop integral partial fractioning is insufficient to eliminate the numerator. This can 
only occur when the reduced graph has a four-point vertex. This is why partial fractioning 
was enough to eliminate the numerator for the box diagram. 
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Eq. (|B23|) at one- loop as described in Ref. [^. The final result for the ^5*1 channel in the 
PDS scheme is 
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-m^ + 2ip^ 

The terms with powers of fi are from a combination of tree, one, and two loop PDS coun- 
terterm diagrams. 

In the ^Di channel the order Q diagrams are shown in Fig. |^. The method used to 
evaluate the box diagram is the same as the ^Si channel, and the only difficult scalar 
integrals that appear are those in Eqs. ( [B16D and (|B20| ). We find 
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This expression agrees numerically with the result in Ref. ||TT|. Fig. yc can be evaluated 
using the result for the NLO one pion exchange ^Si — ^Di diagram 
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APPENDIX C: ORDER Q RADIATION PION CONTRIBUTIONS 



For interactions involving two nucleons it is useful to divide pions into three classes: 
potential, radiation, and soft. This division is analogous to the potential, soft, and ul- 



trasoft regimes |57] devised for calculating non-relativistic diagrams with massless photons 
(NRQED) or gluons (NRQCD) |^8|. To see how the different types of pion arise consider 
evaluating the energy integrals for non-relativistic loop diagrams using contour integration. 
When only residues of nucleon poles are taken, the pions in the graph are potential pions. 
When the residue of a pion pole is taken, the pion is either radiation or soft. A soft pion 
has a momentum which is similar in size to the momentum of the nucleons with which 
it is interacting. A radiation pion exchanges energy with nucleons but does not transfer 
three momentum. Instead, its momentum exchange is governed by a multipole expansion in 



powers of = \Jm.,^/M. Radiation pions are the only type which occur as external particles. 

Loops with only potential or soft pions give functions of p/m^ where p is a nucleon 
momentum. These graphs have a natural power counting in powers of Q ~ p ~ By 
natural power counting we mean that the graph scales homogeneously with Q. On the 



other hand, graphs with radiation pions give functions of p/Qr where Qr = ^jMm^ is the 
momentum threshold for pion production. These graphs have a natural power counting in 
powers of Qr at the scale p ^ Qr This can be seen at the level of the Lagrangian. 



In order to avoid double counting it also necessary to take p ^ Qr when calculating soft 
contributions^. For nucleons with p ^ Qr the three classes of pion are characterized by 
different energy (go) and momentum (g): 

potential go ~ g^/M ~ 

radiation ^Zo ~ 1^1 ~ ""^tt (CI) 



soft go ~ |g I ~ Qr = \J MmT,. 

To implement the KSW expansion, which assumes p ~ m^r, we must expand the result of 
a radiation pion graph in powers of Q. It turns out that the leading Q contribution of a 
radiation pion graph is not determined by the substitution Qr Q^^"^- Instead we will show 



^At p ^ Qr the potential and soft pion propagators should be expanded in TJiT^/Qr. At 
p ~ there may then be factors of vtlt^Ip that must be resummed. See Ref. [^] for an 
explicit example. 



45 



that some radiation pion graphs are enhanced by a factor oil/Q so that an order graph 
can give an order Q^f^'^ contribution. This means that at NNLO the and radiation 
pion graphs need to be considered. In this Appendix we begin by reviewing the power 
counting for pions. The order radiation pion calculation is summarized. We then 
explain how to determine which radiation pion graphs may give an order Q contribution. 
Finally, the order radiation pion graphs which contribute to nucleon-nucleon scattering 
are examined and their order Q contribution is evaluated. 



1. Power counting review 

For A^A^ scattering at NLO the relevant terms in the action are 

S = jdtd?x N^(idt + + Tx\d^t - V' - ml)Tx + -^{N^ a'r^ N){V'7i^) 

-CoOo-D^mlOo + ^O^, (C2) 

o 

where (9o,2 are the four nucleon operators given in Eq. (^. (In this section spin and isospin 
dependence is suppressed since it is not relevant for the rescaling arguements.) To make the 
power counting in this action manifest it is useful to rescale the coordinates and fields in a 
manner similar to the rescaling done in NRQCD p7|J^5| -pT| . The power counting is facilitated 
because factors of p = Mv, and M are made explicit. For the nucleon-pion Lagrangian parts 
of this rescaling were carried out in Ref. and further discussed in Ref. We begin 



by rescaling the coordinates in a manner appropriate to the potential regime and rescaling 
the fields to keep the kinetic terms invariant: 

X T 



X 



Mv Mv 



t=TT^^ ^(^' ^) = (^^)'''' , vr(x, t) = Mv^/' 7rp(X, T) . (C3) 



The coefficients of four nucleon operators will also be rescaled to take into account the KSW 
power counting which is appropriate for large S-wave scattering lengths. Using the PDS O 



or the OS scheme 25,21] and taking fi = Mv gives 



where Cq, C2, -D2 are order v^. This gives the following rescaled action for the potential 
regime 
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Sp = jdTd^X Nl(idT- 
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+ (47r) Mv I ^2(^)' [iVp^iVp]2 + C^iNlNlNpVx + /i.e. 



(C5) 



+ V47r 



Eq. ( |U5D reproduces some familiar features of the power counting. In the nucleon kinetic 
term the Bt and Vj terms are the same order. In the potential pion kinetic term the term 
is down by f ^ and is therefore treated perturbatively. Furthermore, the Vj- and terms 
are the same size for v = itLt^/M ~ 0.15 Thus, p ~ is the natural power counting 
scale when calculating graphs with only potential pions. The Cq interaction term is the same 
size as the nucleon kinetic terms and therefore must be treated non-perturbatively. Each 
potential loop gives a factor of 1/ (47r) which will cancel against factors of (47r) multiplying 
interactions terms like Cq. Insertions of C2 or D2 are suppressed by Mv/K = 771^^/ A ~ 1/2 
and are therefore treated perturbatively. Finally, we see that the exchange of a potential pion 
involves the insertion of two A^A^vr vertices and is suppressed by Mv/Aiy]\f = 111^,/ K^n = 0.47 
where A^^ = {S^p)/{Mg\) ~ 300 MeV. 



In the radiation regime the time coordinate has the same scaling as in Eq. 
spatial coordinate has a different rescaling. 



X 



Tr{x,t) = Mv^ Trr{Xr,T) . 



but the 



(C6) 



The rescaled radiation pion kinetic term is then 
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SXr dT Txl 
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For radiation pions the derivative terms are the same size as the mass term for a different 



value of f, namely Vr = \Jm.,^/M. For v = v^. the radiation pion energy and momentum 
are order m^r. This v corresponds to nucleon momenta p ^ Qr = \fMrnl which is the pion 
production threshold. At these momenta the power counting for graphs with radiation pions 
is straightforward When performing calculations at these momenta the terms in the 

action should be scaled up0 to /i ~ Mv^- The XX-n^ interaction term is 



^We ignore the running of the physical XXtt coupling QAifJ') because its ln(/i) dependence 
is down by Q^. 
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V2 47r/ 

Since the nucleoli and radiation pion fields have a different spatial coordinate we must 



perform a multipole expansion to make the Vr counting manifest, 



(C9) 



Therefore, a nucleon emitting a radiation pion will not have its three momentum changed. 
From Eq. ( |C5| ) we see that each radiation pion vertex comes withQa factor of Mf^/(47r/) = 
m^r/A^. For evaluating radiation pion graphs we take p ~ ~ Qr and have the following 
power counting rules: 

radiation pion propagator 
nucleon propagator 
axial pion — nucleon coupling 
radiation measure 
potential measure 



M/Ql 
Ql/M 
d'q ~ Ql/M' 



(CIO) 



At momenta of order m,r, the mass term in Eq. 
kinetic {d'r 



is enhanced by relative to the 



V|^,) term. For p <^ Qr me see that radiation pions could be integrated out 
in a similar fashion to integrating out W bosons for momenta p -C Mw- Matching onto a 
low energy theory would absorb radiation contributions into local operators. However, this 
will not be done since the matching gives ttIt^ dependence to the coefficients of four nucleon 
operators, yielding a low energy theory without a chiral power counting. Instead, radiation 
pion graphs will be expanded in /Ql Q/M for p ~ m^), and only the order Q piece of 
the radiation pion graphs will be included in our calculation. 

Finally, consider the soft regime Here the spatial coordinate has the same scaling 
as in Eq. (|C3|), but the time coordinate has a different rescaling 
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Mv 



7r{x,t) = (Mv) iTs{X,Ts). 
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The soft pion action is 

= I d^X dT, tt! 
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^Note that since each radiation loop gives a factor of l/(47r)^ we have pulled a (47r) out 
front in the NNtt^ vertex in Eq. 
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where Ns{Ts,X) = Np{T,X). With this rescahng the nucleon action is 

J d^XdT, Nl{idT, - v^)Ns . (C13) 

Therefore, when a nucleon appears in a soft loop the kinetic energy term is treated pertur- 
batively making the propagator static. From Eq. (|C12| ) we see that the power counting of 
soft loops is simplest for v ~ itlt^/M or p ~ m,r- Unfortunately, this makes the soft pion 
modes appear at the same energy and momentum as the radiation pion modes (i.e. ~ ^,7^). 
Therefore, calculating with radiation pions at p ~ Qr and soft pions at p ~ m,r may result 
in double counting. This problem can be avoided by using v = Vr for both radiation and 
soft pions and then scaling down to ~ itlt^/M. An explicit example of this procedure is 
worked out in Ref. [^. Examples of soft diagrams are shown in Fig. p. These diagrams are 
order (even when dressed with Co bubbles) and therefore will not be discussed further. 



2. The order Q part of the order radiation pion graphs 



The order Q^/{M^K^ radiation pion graphs shown in Fig. were calculated^ in 
Ref. [p9|. It is instructive to look at the result of evaluating some of these diagrams: 
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r?V\3/2^ / E 



where fl"^ = 7ie~'^'^ fi"^, and Ji and I2 are hypergeometric functions given in Ref. |29|. The 
1/e poles are cancelled by insertions of a D2rn^ counterterm. The leading order amplitude 
A-i ~ l/(Mp), so we see that Eq. ( |C14|) has terms proportional to 



~M 



Mp 



and 



5/2 



M^/^p'^ 



(C15) 



^The graphs in Fig. |T3a,b and the field renormalization are affected by performing the spin 
and isospin traces in n dimensions, so a) and b) in Eq. (|C14|) differ from Ref. |2^. However, 
the sum of graphs in Eq. ( P16|) is unaffected. 
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a) b) 



f) 5M 5Z 

FIG. 13. Leading order radiation pion graphs for A^A^ scattering. The wavy Hnes are radiation 
pions and 5M, 6Z are the mass and field renormahzation counterterms. There is a further field 



renormalization contribution that is included in Eq. (C16). 



For p ^ Qr these terms scale as Qf/M^, as anticipated by the power counting. At p ~ 
iTT'TT ~ Q, these terms scale like {Q/M)^/'^,Q/M, and {Q/MY^'^ respectively. The graphs 
which give rise to the Q^^^{Q) corrections have two (one) external bubble sums. By external 
bubble sums we mean bubble sums that do not appear inside radiation loops. External 
bubble sums go like 1/p, which scales like 1/Qr at p ~ but 1/Q at p ~ m^. So for each 
external bubble sum, the graph picks up an additional Qr/Q ~ <5^^^^ upon scaling from 
p ~ Qr to p ~ m^r. Terms which scale like Q^^^ at p ~ are actually larger than NNLO 
in the Q counting. The Q^^"^ contributions come from graphs b), e) and /), and cancel when 
these graphs are added together. 

In the ^5*0 channel the sum of all graphs in Fig. Il3 is ||29||: 



- + ln At 



3 ^m. 



where the ln(/i) dependence in Eq. (|C16| ) is cancelled by a ln(/i) in D2 ° (/i). The sum 
of the diagrams turns out to be much smaller than anticipated by the power counting. 
For p ~ Qr, the first term is suppressed by a factor of ~ l/Qrll/a^'^^'^^ — l/a*-^"^^^], the 
second by ~ 1 / Ql[l / a^' ^''^ - l/a('^^^]\ This suppression occurs because the radiation pions 
couple to a charge of Wigner's SU{4) symmetry ||2^, which is a symmetry of the leading 
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order Lagrangian in the limit a*^^'^°\ a^^'^^^ 



oo or a 



The order Ql 



radiation pion graphs are therefore a small correction to the S-wave scattering amplitude. 
Furthermore, to order Q the Qf graphs simply give an additional contribution to the (3 
constant that appears in Eqs. (p!6D, 
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The /i dependence in C^^^ is cancelled by /i dependence in D2 



The result in the ^5*1 
channel is obtained from Eqs. ( |(J16D and ( |(J17D by switching the ^Sq and ^5*1 labels. 



3. Scaling radiation contributions from Qr to m„ 



Since we are interested in the power counting for p ~ it is important to know how 
big a radiation pion graph may get when p is lowered from Qr to m^r. The graphs have 
pieces that scale as Q^^'^,Q, Q^^'^, . . ., for p ~ as discussed in the previous section. In 
order to know which radiation pion graphs to include at a given order in the KSW power 
counting, we must know the size of the leading term in the Q expansion of a graph for 
p ~ m^r. In this section we will prove that an order Q^. calculation is sufficient to determine 
the order Q^l"^'^ result. 

To see this first consider the Q expansion of pcot5 in the ^5*0 channel: 
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pcot S is real and an analytic function of p"^ near p = 0. This will be true order by order in 
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where the fk are real functions of p which are analytic about — 0. The general form of 
a higher order amplitude is powers of A-i multiplied by functions of p. The crucial point 
is that the function multiplying the [^-i]^ is the only new contribution. The coefficient of 
[.4_i]"*, m > 2, is determined by lower order amplitudes. The graphs giving the m > 2 
contributions are "Co reducible" by which we mean that they fall apart when cut at an A-i 
vertex. 

This generalizes to the Qr expansion of radiation pion graphs, the only difference being 
that the radiation pion contribution starts out at Qf, while the potential pion starts out at 
(5°. ^ Qr radiation pion correction to the amplitude will be of the form: 

Ak = [A-,]'fk,2 + [A-,ffk,3 + ... + [A-,]'-'fk,k-i. (C20) 

Again, the fk,m are real and analytic about p^ = and all the fk,m except for 2 will be 
determined from lower order amplitudes. Since Ak ~ Qr and A-i ~ l/{Mp), fk,2 ~ Qr"*"^ 
for p Qr- To understand how fk^2 scales with Q as p is lowered to 171^^, note that without 
loss of generality, fk,2 can be written as 

where the ellipses denote momentum dependence that involves scales other than Q^, and 
A — A^, A, or M. For p ~ m,r the ellipse denote dependence on the dimensionless variables 
p/rriT^^ ap, and p/ A. For p ~ m^, p/y/Mrn^ ~ {Q/My/'^ and the function fk,2 can be 
expanded in its first argument: 



^-l/fc,2 - ^-1 k-2 Jk,2 lU, . . .) 

A2 A 



1/2- 

1 + " ' 



Q_V 

m) 



(C22) 



Therefore, the new contribution at Q^ scales hke Q^l'^~^ (plus subleading terms) forp ~ 7x1^^. 
This is consistent with the result of the Q\ calculation, where the largest contributions from 
individual graphs scaled as Q^l"^ . A cancellation between graphs resulted in this contribution 
vanishing. The remaining terms scale as Q, Q^^^, .... 

Next we consider contributions to the amplitude from Cq reducible graphs. If a Cq 
reducible graph is obtained by joining j Cq irreducible graphs where the j'th graph scales 
as Q'^i at p ~ m-n-, then the Cq reducible graph scales as 

g^-i+E-=i"i. (C23) 
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x^GXD*^ x^ooO*^ x*0*O*^ 



FIG. 14. Example of order graphs that have three external bubble sums. 

For example, the order graphs in Fig. |T^ are each obtained by joining a potential 
pion graph with a radiation pion graph. The radiation graphs scale as Q^^"^ for p ~ 
so the individual graphs in Fig. [14] scale as Q^^^"^ for p ~ m^r. No Cq irreducible graphs give 
an order Q contribution, 
bmce 

_ g ^ _ 4^ Q4 radiation pion graphs can have a contribution that is 
NNLO for p ~ m^. This calculation is taken up in the next section. Note that a calculation 
of the order graphs would be necessary to determine the order Q^/^ terms. 



4. The order Q part of the order radiation pion graphs 

The order radiation pion contributions come from graphs that have one radiation 
pion, an arbitrary number of Cq's, and one insertion of a C2P^, D2'm'^, or G2 operator or 
one potential pion. The coefficient G2 multiplies a four-nucleon operator that couples to the 
axial pion current, 

C='-G2 [iV^P/^^iV]t [iV^p/^V,(e5,-e^ - ed,i)N] + h.c. . (C24) 

Note that due to the hermitian conjugate this operator is the same for s = ^Sq and s = ^Si. 
Power counting these graphs gives Q^/(M^A^A), i.e. they are suppressed by Qr/A relative to 
the leading radiation pion graphs in Fig. Note that Qr = 360 MeV, so for A < 360 MeV, 
the Qr/-^ expansion does not converge. If this is the case then the radiation pion contribution 
is not calculable. This is true of radiation contributions even when we scale down to p ~ m^r. 
To make the radiation contributions calculable we must have a power counting for the pure 
potential contributions that works for p ~ Qr- One possible resolution is to ignore radiation 
pion contributions since at low momenta the radiation pions can be integrated out. However, 
this makes the coefficients of four nucleon operators depend on 171.,^ in a non-trivial way. We 
will proceed by computing the radiation contribution which is formally order Q in both 
S-wave channels even though the size of the spin triplet potential diagrams in section [111 B 
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indicate that a modification of tlie power counting is likely necessary to obtain a convergent 
expansion in this channel. Since the C2P^ and D2Tn^ operators and potential pion exchange 
do not respect Wigner symmetry, there will be no suppression by factors of l/{aQr) at this 
order. 

For calculational purposes it is useful to define offshell amplitudes for S-wave transitions, 
N{pi)N{p2) —>■ N{p3)N{p4), induced by 4-nucleon operators. These amplitudes are equal to 
a sum of Feynman diagrams where the equations of motion have not been used. They can 
be treated as vertices and inserted inside loop graphs, which greatly reduces the number of 
order Qj diagrams. The offshell order 1/Q amplitude in n dimensions is 



where 



471" r(l - 72/2) //ix3-n 

and E is the center of mass energy 

E-E, + E2 -E, + E,- 4^ . (C27) 

At order the NN amplitude has contributions from the four nucleon operators C2, -D2, 
and Cg''''. In n dimensions the offshell amplitude for C2 graphs is 

- 9. -XX + XX -XXX 



Cn ^ ^ Cn V 8 J 



Note that C2(/i)/Co(/^) and therefore the offshell C2 amplitude are /i dependent. The onshell 
amplitude is /i independent using the order l/Q^ part of the beta function in Eq. (^) since 
the term proportional to C2/C0 vanishes by the equations of motion. The offshell amplitude 
that includes the graphs with D2 or Cq^"* vertices is 

£'2X0 ^2,^0 ^2,^0 i^2,CQ 



" (Co) 



] (C29) 
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a) b) 



2>®<C 



c) d) e) 



8^^®tX^ 4)><X>< 



f) . g) . h) 



FIG. 15. Order Qf radiation pion graphs with insertions of C2 and D2. The soHd lines are 
nucleons and the wavy hnes are radiation pions. The black circle denotes the Co bubble sum, 
iA^~^\ and the hatched circle denotes an insertion of iAc2 ^-^D2 given in Eqs. ( P28D and ( |C29D . 



The order Qf radiation pion graphs with insertions of iAc2 ^-^^2 are shown in Fig. [15. 
We find that graphs with insertions of iAD2 gi'^^ contributions that are order or higher. 
The graphs which have an iAoz external to the radiation loop (Fig.|l5|b,f,h) have no order 
Q contribution because the same cancellation that occurs in Fig. Il^,e,f occurs here. Of 



the remaining graphs only Fig. |T5|e with two external bubble sums can have an order Q 
contribution, however inside the radiation pion loop iAD2 ~ ~ Q ('^ot Q^) so this 

graph is order Q^. 



With insertions of iAc2 5 ^^e only graphs in Fig. |T5| which do not give an order Q con- 
tribution are a) and c). Diagrams d) and g) give a non-zero order Q contribution even 
though they have only one external bubble sum. The order Q contribution comes from the 
/i dependent part of iAc2- Since these graphs have one external bubble sum they are 
expected to be ~ Q^^"^. However, with {pi — ~ Ql or ME ~ Q^, and jj, Q the 11 
dependent part of iAc2 is order Qr/ /i ~ 1/ -/Q. This extra factor^ makes these graphs order 
Q. For the ^5*0 channel, the sum of the order Q contributions from the C2 radiation graphs 
in Fig. O is 



^Note that enhancements by factors of /x do not effect the proof in section C3 since the 
amplitude at a given order in Qr is /z independent. 
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FIG. 16. Order Qf. radiation pion graphs with insertions of G2. 



• Arad 



P 



2 Mml (m^ 



47r 



C. 



(C30) 







The corresponding amplitude in the ^Si channel is obtained by exchanging the ^5*0 and ^5*1 
labels in Eq. ([USOl) . 

The result in Eq. ( |C30|) is /i dependent. For the term proportional to the /z de- 



pendence is cancelled by a radiation contribution to the beta function for D2{fi). The /i 
dependence of the term is cancelled by 6*2 (z^)- To calculate the PDS beta function for 
G2(/i) we consider graphs with G2 dressed with Cq bubbles on one side and Cq bubbles 
on the other. We also consider graphs with a C2 vertex next to a NNn vertex: 

C2 

(C31) 

When the derivatives in O2 act on the nucleons on the right there is a piece in which the 
numerator cancels the propagator exactly. This piece has the same form as a NNNNir 
vertex and contributes to the beta function for G2 when dressed with Cq bubbles. We find 




A 



^(i5o)^(35i) 



which has the solution 



G2W = K,G 



Ma 



(C32) 



(C33) 



47r An 

with Kg the constant of integration. Eq. ( |(J33| ) gives G2{fi) ~ V/^^- Because G2 is the 
coefficient of the four nucleon coupling to the axial current it has the same renormalization 
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8 



4 



-/V 



FIG. 17. Order Qf radiation pion graphs with one potential pion. 



group equation as the weak axial four nucleon operator[3 considered in Ref. [|T3|. Using the 
scahng of G2 to power count the diagrams in Fig. 16 we find that they are order Qf radiation 
pion graphs. For the ^5*0 channel, the sum of the order Q part of the diagrams in Fig. |16 



IS 



• AT ad 



G2 



An 



(C34) 



The result in the ^Si channel is obtained by interchanging the labels ^5*0 and ^Si. Using 
Eq. ( |C33| ) we see that the /x dependence of the term in Eq. ( |C30| ) is cancelled by the 
term in Eq. (|U34|) . The /i dependence of the term in Eq. (|(J34|) is again cancelled by 



DM. 

The final order Qf. diagrams that we must consider are those with one potential pion, 



one radiation pion and an arbitrary number of Cq's shown in Fig. |T^. Here we find that only 
graphs with two external bubble sums can give an order Q contribution. Graphs in Fig. |l^ 
with three external bubbles sums are Co reducible and do not give order Q contributions. 
Furthermore, of all the diagrams in Fig. ^ with two external bubble sums only the three 



shown in Fig. |l^ give an order Q contribution. These are the graphs in which the potential 
pion exchange is inside the radiation pion loop. This ensures that all potential loop momenta 
in the graph see the scale Qr which is necessary for the graph to give an order Q contribution. 
So, for example, we find that the graph in Fig. IT^ does not give an order Q contribution. 



The diagrams in Fig. |T8| look somewhat daunting since they involve a three, four, and five 
loop calculation. Nevertheless, their order Q contribution can be evaluated analytically. 



Adding up the the order Q part of the diagrams in Fig. jTBI gives 



lOjn Ref. [H G2 was denoted by Li^a- 
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FIG. 18. The three order Qf. radiation pion graphs with one potential pion that give an order 
Q contribution for p ~ tjItj. 




FIG. 19. An order Qf radiation pion graphs with one potential pion that does not give an order 
Q contribution. 



It is interesting that the order Q part of these graphs is not multiphed by a nontrivial 
function of p/mT^. At one higher order, Q^^"^^ the diagrams in Fig. |l^ will give a result which 
involves a function of p/mT^. 

The results in Eqs. (|C30| , |C34|C35| ) give the complete order Q contribution from order Q'^ 



graphs. It is interesting to note that all contributions are equal to a constant times 
Therefore they simply give an additional contribution to the constant that appears in 
Eqs. ([T|) and 



The result for (^'^^ is the same in the ^Si channel. The n dependence in (^"^^ is cancelled by 
fi dependence in D2. 

The complete order Q contribution from radiation pion graphs is the sum of Eqs. ( |C17|) 
and (IC36D : 



cr' = d'^ + ct^ ■ (C37) 

The order Q radiation pion contribution to the D2 "^"'^ beta function in Eq. (P) is 
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The contribution to the D2 beta function in Eq. (|^) is obtained by switching the spin 
singlet and triplet labels. 

For a consistent radiation pion calculation at momenta p ~ it would be necessary 
to keep all powers of p/Qr in computing the diagrams. If the Qr/A expansion were 
convergent then the magnitude of this radiation contribution would be small; being down 
by Ql relative to the leading order amplitude. For Qr/A > 1 it is necessary to modify the 
power counting for the potential diagrams to increase the scale A before the radiation pion 
power counting will yield a convergent series. Since diagrams with one radiation pion are 
suppressed by m^/A^ relative to any leading order amplitude these contributions are always 
likely to contribute at the few percent level. 



59 



REFERENCES 



[1] J.F. Donoghue, E. Golowich and B.R. Holstein, Dynamics of the Standard Model (Cam- 
bridge University Press,, 1992). 

[2] E. Jenkins and A. Manohar, UCSD-PTH-91-30; E. Jenkins and A. Manohar, Phys. Lett. 
B255 (1991) 558, ihidB 259 (1991) 353. 

[3] V. Bernard, N. Kaiser, and U. Meissner, Phys. Rev. Lett. 67 (1991) 1515; Phys. Lett. 
B 319 (1993) 269; Int. J. Mod. Phys. E4 (1995) 193. 

[4] S. Weinberg, Phys. Lett. B251 (1990) 288; Nucl. Phys. B363 (1991) 3; 

[5] C. Ordonez and U. van Kolck, Phys. Lett. B291 (1992) 459; C. Ordonez, L. Ray and 
U. van Kolck, Phys. Rev. Lett. 72 (1994) 1982; U. van Kolck, Phys. Rev. C49 (1994) 
2932; G.P. Lepage, |nucl-th/9706"029| ; T-S. Park, K. Kubodera, D-P. Min and M. Rho, 
Nucl. Phys. A646 (1999) 83; D.R. Phillips and T.D. Cohen, ^ucl-th/9906WTi See the 
review, U. van Kolck, |nucl-th /99020 15| and references therein. 

[6] C. Ordonez, L. Ray, and U. van Kolck, Phys. Rev. C53, (1996) 2086, 

[7] E. Epelbaum, W. Glockle, and U.G. Meissner, |nucl-th/9910064| , 

[8] D. B. Kaplan, M. J. Savage, and M. B. Wise, Phys. Lett. B424 (1998) 390, 

[9] D. B. Kaplan, M. J. Savage, and M. B. Wise, Nucl. Phys. B534 (1998) 329, 
[10] U. van Kolck, |hep-ph/971122^ , 

[11] N. Kaiser, R. Brockmann, and W. Weise, Nucl. Phys. A625 (1997) 758. 
[12] E. Epelbaum and U.-G. Meissner, ^iucl-th/9903046t X. Kong and F. Ravndal, |ieF 
I ph/ 99035231 ; Phys. Lett. B450, (1999) 320; X. Kong and F. Ravndal, |nucl-th/ 990206^ ; 
|nucl-th/9904"06^ ; D. B. Kaplan, M. J. Savage, and M. B. Wise, Phys. Rev. C59 (1999) 
617; J.-W. Chen et.al, Nucl. Phys. A644 (1998) 221; M. J. Savage, K. A. Scaldeferri, 
and M. B. Wise, |nucl-th/98lT029| ; J.-W. Chen, G. Rupak, and M. J. Savage, [riucl^ 



th/99050U^ ; J.-W. Chen, H. W. Griesshammer, M. J. Savage, and R. P. Springer, Nucl. 
Phys. A644, 245 (1998); J.-W. Chen, |nucl-th/981002Tt D. B. Kaplan, M. J. Savage, 
R. P. Springer, and M. B. Wise, Phys. Lett. B449 (1999) 1; M. J. Savage and R. P. 
Springer, Nucl. Phys. A644 (1998) 235; P.F. Bedaque and H.W. Griesshammer, [nucl 
th/99070771 . 



[13 
[14 

[15 
[16 
[17 
[18 
[19 

[2o; 

[21 
[22 
[23 
[24^ 
[25 
[26^ 



Malcolm Butler and Jiunn-Wei Chen, |nucl-th / 9905059] , 

M. J. Savage, |nucl-th / 9905009| , Proceedings of the INT Workshop on Nuclear Physics 
with Effective Field Theory, Seattle, WA 25-26 Feb 1999. 
M. Binger, |nucl-th/99010T2| , 
G. Rupak and N. Shoresh, |nucl-th/ 9902077 , 
G. Rupak and N. Shoresh, |nucl-th/9906W7 , 



T. Mehen and I. W. Stewart, |nucl-th/99060T0 



S. Fleming, T. Mehen amd I. Stewart, |nucl-th/9906056 



V.G.J. Stoks, et.al., Phys. Rev. C48 (1993) 792; V.G.J. Stoks et.al, Phys. Rev. C49 

(1994) 2950, |nucl-th/9406039| . (cf. |http://nn-online.sci.kun.nl/NN/|) 

T. Mehen and I. W. Stewart, Phys. Rev. C59 (1999) 2365, 

E. Wigner, Phys. Rev. 51, 106, 947 (1937). 

T. Mehen, I. W. Stewart, and M. B. Wise, |hep-ph/ 9902370 , 

J. Gegelia, ^ucl-th/ 980205^ . 

T. Mehen and I. W. Stewart, Phys. Lett. B445, (1999) 378, 



M.C. Birse, J.A. McGovern, and K.G. Richardson, |hep-ph/ 9807302 



60 



[27' 

[28^ 
[29 
[30 
[31 
[32 
[33^ 
[34 
[35^ 
[36^ 
[37^ 
[38^ 



I. W. Stewart, Ph.D. thesis, |hep-ph/9907448 



J.-W. Chen, G. Rupak, and M. J. Savage, Nucl.Phys. A653 (1999) 386, 
T. Mehen and I. W. Stewart, |nucl-th/990106l| , 
J. Steele and D.B. Kaplan, |nucl-th/9905027| . 
T.D. Cohen and J.M. Hansen, Phys. Rev. C59 (1999) 13; |nucl-th/9901065 
D.B. Kaplan, private communication, 

H.P. Stapp, T.J. Ypsilantis, and N. Metropohs, Phys. Rev. 105 (1957) 302, 
N. Kaiser, S. Gerstendorfer, and W. Weise, Nucl. Phys. A637 (1998) 395, 
J.W. Chen, G. Rupak, and M. J. Savage, |nucl-th/ 990500^ , 
D.R. Phillips, G. Rupak, and M.J. Savage, |nucl-th/990805l , 



M. Luke and A. Manohar, Phys. Rev. D55 (1997) 4129, 

P.P. Bedaque, H.W. Hammer and U. van Kolck, Nucl. Phys. A646, (1999) 444; P.P. Be- 
daque, H.W. Hammer and U. van Kolck, |nucl-th / 9906032] . P.P. Bedaque and U. van 
Kolck, Phys. Lett. B428, (1998) 221; P.P. Bedaque, H.W. Hammer and U. van Kolck, 
Phys. Rev. C58, (1998) R641. 
[39] P.P. Bedaque and H.W. Griesshammer, |nucl-th/9907"077| . 

[40] F.V. Tkachov, Phys. Lett. BlOO (1981) 65; K.G. Chetyrkin and F.V. Tkachov, Nucl. 

Phys. B192 (1981) 159, 
[41] G. Passarino and M. Vehman, Nucl. Phys. B160 (1979) 151; G. Weiglein, R. Scharf, 

and M. Bohm, Nucl. Phys. B416 (1994) 606, 
[42] O.V. Tarasov, Phys. Rev. D54 (1996) 6479; Nucl. Phys. B502 (1997) 455. 
[43] R. Mertig and R. Scharf, Comput. Phys. Commun. Ill (1998) 265. 
[44] R. Harlander and M. Steinhauser, |hep-ph/9812357| , 
[45] E. Braaten and A. Nieto, Phys. Rev. D51, (1995) 6990. 
[46] A.K. Rajantie, Nucl. Phys. B480 (1996) 729. 
[47] M. Beneke and V.A. Smirnov, Nucl. Phys. B522 (1998) 321, 

[48] W.E. Caswell and G.P. Lepage, Phys. Lett. B167 (1986) 437; G.T. Bodwin, E. Braaten, 

and G.P. Lepage, Phys. Rev. D55 (1997) 1125. 
[49] B. Grinstein and I. Rothstein, Phys. Rev. D 57 (1998) 78; P. Labelle, Phys. Rev. D58 

(1998) 093013; M. Luke, A.V. Manohar, and I.Z. Rothstein, |hep-ph/9910209| , 
[50] M. Luke and M.J. Savage, Phys. Rev. D57 (1998) 413; 
[51] H.W. Griesshammer, Phys. Rev. D58 (1998) 094027 and |hep-ph/9810235| , 
[52] H.W. Griesshammer, |hep-ph/ 98042511 



61 



